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Abstract 

We consider the two-nucleon system at next-to-next-to-next-to-leading order (N 3 LO) in chiral 
effective field theory. The two-nucleon potential at N 3 LO consists of one-, two- and three-pion 
exchanges and a set of contact interactions with zero, two and four derivatives. In addition, one 
has to take into account various isospin-breaking and relativistic corrections. We employ spectral 
function regularization for the multi-pion exchanges. Within this framework, it is shown that the 
three-pion exchange contribution is negligibly small. The low-energy constants (LECs) related 
to pion-nucleon vertices are taken consistently from studies of pion-nucleon scattering in chiral 
perturbation theory. The total of 26 four-nucleon LECs has been determined by a combined fit to 
some np and pp phase shifts from the Nijmegen analysis together with the nn scattering length. 
The description of nucleon-nucleon scattering and the deuteron observables at N 3 LO is improved 
compared to the one at NLO and NNLO. The theoretical uncertainties in observables are estimated 
based on the variation of the cut-offs in the spectral function representation of the potential and 
in the regulator utilized in the Lippmann-Schwinger equation. 
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1 Introduction 



Since the seminal work of Weinberg [1] to derive the forces between two, three, . . . nucleons from 
chiral effective field theory, there has been a flurry of activities to work out the consequences of such 
an approach, to improve that scheme or to construct alternatives, for reviews see [2, 3]. Here, we 
will be dealing with a modified Weinberg scheme, in which pions are treated nonperturbatively and 
the power counting is applied to the nucleon-nucleon potential. The potential consists of one-, two-, 
. . . pion-exchanges (1PE, 2PE, . . .) and a string of contact interactions with an increasing number 
of derivatives (zero, two, four, . . .) that parameterize the shorter ranged components of the nuclear 
force (the precise framework is specified in more detail below). Such an approach has a variety of 
advantages over more conventional schemes or phenomenological models. First, it offers a systematic 
method to improve calculations by going to ever increasing orders in the power counting and it 
allows to give theoretical uncertainties. Second, one can consistently derive two- and three-nucleon 
forces (see e.g. [4]), which has never been achieved before and paves e.g. the way for a new look 
at the problem of nuclear matter. Third, nucleon and nuclear properties can be calculated from one 
effective Lagrangian, which is of particular importance if one intends to extract neutron properties 
from (electromagnetically induced) measurements on light nuclei in a controlled theoretical way. In 
this paper, we present the nucleon-nucleon potential at next-to-next-to-next-to-leading order (N 3 LO) 
in the chiral expansion, extending our earlier work, and we apply this potential to observables in two- 
nucleon systems. Our work differs from the one of Entem and Machleidt (EM), who first presented an 
N 3 LO potential in Ref. [5], in various ways. First, they use the two-pion exchange contributions based 
on dimensional regularization, which have a very singular short-range behavior. We employ spectral 
function regularization, which allows for a better separation between the long- and short-distance 
contributions. Second, EM present results only for one choice of the cut-off necessary to regulate the 
high-momentum components in the Lippman-Schwinger equation to generate the scattering and the 
bound states.* 5 We perform systematic variations of this cut-off and other parameters which allows 
us to give not only central values but also theoretical uncertainties. Third, our treatment of the isospin 
breaking effects differs from the one of EM (which is based on our earlier work [7]). Fourth, we employ 
a relativistic version of the Schrodinger equation, which allows to calculate consistently relativistic 
corrections also in three and four nucleon systems. Other less significant differences will be discussed 
in due course. We believe that with the material presented here an important step has been made 
to put precision calculations in nuclear physics on a firm theoretical basis, which not only allows to 
readdress many issues that have already been investigated in quite a detail but will also open new 
areas of testing chiral dynamics in few-nucleon systems or shed more light on the issue of the nuclear 
forces in the limit of vanishing quark masses (see [8, 9, 10, II] for earlier work on that topic). 

Our manuscript is organized as follows. In sec. 2, we explicitly give the potential at N 3 LO. The 
contributions up to NNLO have already been extensively discussed in [12, 13, 14, 15, 16, 17]. The 
N 3 LO corrections due to two- and three-pion exchange have been derived recently in [18, 19, 20] 
using dimensional regularization to regularize divergent loop integrals. Throughout this work, we use 
a different regularization scheme for the potential, namely the spectral function regularization (SFR). 
This approach has been recently proposed and successfully applied at NLO and NNLO [17, 21]. As 
demonstrated in these references, the SFR scheme allows to significantly improve the convergence of 
chiral Effective Field Theory (EFT) for the two-nucleon system. We also give an overview of various 
isospin-breaking interactions including electromagnetic forces and discuss the regularization procedure 
necessary to render the (iterated) potential finite. In section 3 we deal with the scattering equation. 
In order to account for the relativistic corrections to the nucleon kinetic energy, we have decided to 

#5 The results for two different choices of the cut-off in the Lippmann-Schwinger equation are shown in [6]. 
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use the Lippmann-Schwinger equation with the relativistic expression for the kinetic energy. Such an 
approach can naturally be extended to few-nucleon systems and to processes with external probes. 
We also discuss how to cast the relativistic Lippmann-Schwinger equation into a nonrelativistic form, 
which might be useful in certain applications. Various deuteron properties are considered in section 4 
using both relativistic and nonrelativistic Schrodinger equations. The fitting procedure to determine 
the low-energy constants (LECs) and the accuracy of the fits are detailed in sec. 5. Results for 
phase shifts and the deuteron (bound state) properties are displayed and discussed in sec. 6. Our 
findings are summarized in sec. 7. The appendices contain details on the kinematics, the partial wave 
decomposition, the momentum space treatment of the Coulomb interaction and the effective range 
expansion. 



2 The two— nucleon potential at N 3 LO 
2.1 General remarks 

Before going into details of calculations, we would like to make certain general remarks. As already 
pointed out in the introduction, we strictly follow the scheme suggested by Weinberg [1]. In this 
approach one uses the EFT technique to derive nuclear forces from the most general (approximately) 
chiral invariant effective Lagrangian. The NN S-matrix is obtained via (non-perturbative) solution of 
the Lippmann-Schwinger (LS) equation. In most practical calculations (including the present one), 
the later step can only be performed numerically. 

Starting from the most general chiral invariant effective Hamiltonian density for pions and nucleons 
one can derive energy-independent and hermitean nuclear forces by a variety of methods including 
the method of unitary transformation, see e.g. [10]. The resulting nucleonic forces are ordered by the 
power of the generic low-momentum scale Q related to the three-momenta of nucleons, the pion mass 
and typical four-momenta of virtual pions: 



V ~ O 



i_ (Q_ 



(2.1) 



where F n is the pion decay constant and A x is the chiral symmetry breaking scale or, more generally, 
the smaller of the chiral symmetry breaking scale and the scale Alec associated with the LECs 
accompanying four-nucleon contact interactions (as discussed below). The power v for a given diagram 
can be calculated using the rules of dimensional analysis [1] 

v = -2 + 2E n + 2(L - C) + V i A i > ( 2 - 2 ) 

i 

where E n , L, C and Vi are the numbers of nucleons, loops, separately connected pieces and vertices 
of type i, respectively. Further, the quantity Aj, which defines the dimension of a vertex of type i, is 
given by 

At = di + -m - 2 , (2.3) 

with di the number of derivatives or M n insertions and ni the number of nucleon lines at the vertex 
i. One has A, > as a consequence of chiral invariance. This leads to u > for connected diagrams 
with two and more nucleons. One also recognizes that the graphs with loops are suppressed and that 
(n + l)-nucleon forces appear at higher orders than the n-nucleon ones. We note, however, that 
the formula (2.2) does not apply to a specific sort of diagrams, sometimes referred to as reducible, 
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whose contributions are enhanced due to the presence of anomalously small energy denominators 
resulting from purely nucleonic intermediate states. Such reducible diagrams are responsible for the 
non-perturbative aspect in the few-nucleon problem and must be summed up to infinite order. They, 
however, do not contribute to the nuclear potential and result from iteration of the potential in the 
Lippmann-Schwinger equation. 

It remains to specify our way of counting the nucleon mass. In the single-nucleon sector it appears 
to be natural to treat the nucleon mass m in the same way as the chiral symmetry breaking scale 
A x ~ 1 GeV. As argued in [1], in the few-nucleon sector consistency requires that the nucleon mass is 
considered as a much larger scale compared to the chiral symmetry breaking scale. If one adopts the 
counting rule m ~ A x , no nonperturbative resummation of the amplitude is required from the point 
of view of the chiral power counting. In this work we adopt the counting rule m/Q ~ (A x /Q) 2 , which 
has also been used in [12]. 

In the following sections we will discuss various contributions to the NN potential up to N 3 LO including 
isospin-breaking corrections. 

2.2 Contact terms 

We consider first the contact terms of the two-nucleon potential. To the accuracy we are working, 
the potential in the center-of-mass system (cms) for initial and final nucleon momenta p and p ', 
respectively, takes the form:* 6 



eq. (2.2). Notice that the contact operator basis in eq. (2.4) represents just one particular choice among 
many others. One could equally well use another set of 24 independent contact operators including 
for instance terms which contain the product of isospin matrices t\-T2- A one-to-one correspondence 
between different sets of contact operators can be established upon performing antisymmetrization of 
the potential, see ref. [22] for more details. Notice that we have only shown isospin-invariant terms 
in eq. (2.4). Isospin-breaking short-range corrections will be specified below. 

The terms in eq. (2.4) feed into the matrix-elements of the two S-waves ( 1 S'o, 3 Si), the four P-waves 
( 1 Pi, 3 Pi, 3 P2, 3 Po), the four D-waves ( 1 Z) 2 , 3 D 2 , 3 D 3 , 3 Di) and the two lowest transition potentials 

#6 We use the notation of ref. [16] (ref. [6]) for V (0) and V (2) (V (4) ). 
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with q = p ' — pand k = (p+p')/2. The superscripts denote the corresponding chiral order as defined in 
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with p 



and p' = |p'|. The spectroscopic LECs are related to the ones in eq. (2.4) according to 
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D 3P1 = ~{32D 1 -2D 2 + 32D 5 -2D 6 -8D 9 -2D 10 + 48D 11 + 4D 12 -4D 13 -3D 14 ) , 
o 

D 3 P2 = -^-(160Di - WD 2 + 160£> 5 - 10D 6 + 40A) + IOD10 + 16Z>n - 4D 12 + 4L> 13 - D 14 ) , 
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15 
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15 

= D ll = ^( 112D H " 4j °12 " 4D 13 + 7^14 - 16D 15 ) , 
D 2 3D1 - 3S i = D 2 el = ^(16Dn + 4D 12 + 4,0x3 + D 14 ) , 

D 3F2 - 3 P2 = D e2 = -^(16L> n - 4D 12 + 4D 13 - D 14 ) , (2.6) 

These 24 constants are not fixed by chiral symmetry and have to be determined by a fit to data or 
phase shifts and mixing parameters in the corresponding channels. From each of the two S-waves, we 
can determine four parameters, whereas each of the four P-waves and the mixing parameter ei contain 
two free parameters. Further, one free parameter contributes to each of the four D-waves and to the 
mixing parameter e 2 . Of course, we have to account for the channel coupling in the mixed spin-triplet 
partial waves. Once the spectroscopic LECs have been determined, the original Cs, Ct, Ci, . . . , C 7 
and D\ , . . . , D15 are fixed uniquely. 

2.3 One—, two— and three-pion exchange 

Consider now one-, two- and three-pion exchange (3PE) contributions V l7r , V 2lT and V 3lT , respectively. 
At N 3 LO (Q 4 ) in the low-momentum expansion Q they can be written as: 

t/ _ t/(0) , T/ (2) T/ (3) T/ (4) 
V2n ~ V 2n +V 2n +V 2n + • • • , 

V* = V£ ] + .... (2.7) 

Here the superscripts denote the corresponding chiral order and the ellipses refer to Q 5 - and higher 
order terms which are not considered in the present work. Contributions due to exchange of four- and 
more pions are further suppressed: n-pion exchange diagrams start to contribute at the order Q 2n " 2 , 
see e.g. [1]. 

In the following we will give explicit expressions for the individual contributions in eq. (2.7). The 
pion-exchange NN potential in the two-nucleon center-of-mass system (c.m.s) takes the form: 

V = Vc + Ti-t 2 Wc + [Vs + Ti • r 2 Ws] &1-&2 + Wt + Ti • T2 Wt] a± ■ qa 2 ■ q (2.8) 
+ [Vls + ti-t 2 W LS ] i(ai + a 2 ) -(qxk) + [V aL + n ■ t 2 W aL ] a x ■ (q x k)a 2 -(qxk), 
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where the superscripts C, S, T, LS and ah of the scalar functions Vc, W a L refer to central, 
spin-spin, tensor, spin-orbit and quadratic spin-orbit components, respectively. 

The leading order 1PE potential is given by 

Wt {q) - \2fJ q 2 +M 2 ■ (2 - 9) 

At NLO one has to take into account various corrections which result from one-loop diagrams with 
the leading vertices and tree graphs with one insertion of the die- and ^-vertices (in the notation 
of rcf. [23]) from the dimension three Lagrangian C^n and ^4-vertices from the dimension four 
Lagrangian C n . All these graphs lead just to renormalization of the LECs gA, and the pion mass 
M n , see [10] for more details. In addition, one has a contribution from the 1PE graphs with one 
g^-vertex replaced by the dis-vertex from the dimension three ttN Lagrangian. This correction leads 
to the so-called Goldberger-Treiman discrepancy and can be accounted for by the replacement 

9A^9A~ 2d i8 M 2 (2.10) 

in Eq. (2.9). The corrections at NNLO arise from one-loop diagrams with one subleading nNN vertex 
and lead to renormalization of the LEC qa [14]. The corrections to 1PE at N 3 LO are due to two-loop 
diagrams with all vertices of the lowest chiral dimension, one-loop graphs with one subleading vertex 
from C n N or jC n and tree graphs with two subleading vertices or one sub-subleading vertex. After 
performing renormalization of the LECs, one finds the N 3 LO contribution to the 1PE potential to be 
proportional to [24] : 

Wi%) oc ^f/tf) , (2-11) 

where the function F{q 2 ) = aiM% + a 2 M 2 q 2 + a^q 4 can be viewed as the pion-nucleon form-factor. 
The latter does not represent an observable quantity. Expressing the function F(q) as 

F{q 2 ) = (ai -a 2 + a 3 )M* + (q 2 + M 2 )(a 2 M 2 + a 3 (q 2 - M 2 )) (2.12) 

the N 3 LO contribution reduces to a renormalization of the 1PE potential (2.9) and contact interactions. 
Notice that one also encounters an additional correction to the Goldberger-Treiman discrepancy. In 
addition, one has to take into account relativistic l/m 2 -corrections to the static 1PE potential, which 
depend on a particular choice of the unitary transformation in the NN system and the form of the 
scattering (or bound state) equation, see ref. [25] for details. The final expression for the 1PE potential 
adopted in the present work takes the form 

9a\ 2 A p 2 +p' 2 \^ ^Si-qa 2 -q 



VlM = -\M) 1 1 -1S-J ,1 ' S 7W' (2 ' 13) 

where, as in our previous work [21], we take the larger value qa = 1-29 instead of qa = 1-26 in order 
to account for the Goldberger-Treiman discrepancy. This corresponds to the pion-nucleon coupling 
constant g^N = 13.1. 

We now turn to the 2PE contributions. The 2PE potential V$+V$ is discussed in [13, 14, 15, 16] and 
in [12] using an energy-dependent formalism. While dimensional regularization or equivalent schemes 
have been used in [13, 14, 16] to calculate matrix elements of the potential, a finite momentum cut- 
off approach has been applied in [12]. The N 3 LO corrections have been recently obtained by Kaiser 
using dimensional regularization [18]. In the following, we will adopt the SFR method to obtain the 
non-polynomial contributions to the 2PE potential with the short-range components being explicitly 
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excluded, see [17] for more details. The expressions for the 2PE potential in the SFR scheme up to 
NNLO have already been given in [21]. To keep the presentation self-contained, we give below the 
corresponding (non-polynomial) terms at NLO 



= -^I^ LA ^){ 4M ^A- 4^-1) +q 2 (23g\ -10^-1) + ^^ 



and at NNLO: 
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(2.14) 



(2.15) 



where the NLO and NNLO loop functions L A (q) and A A (q) are given by 

to K 2 lo 2 + q 2 s 2 + 2Agws 



L a (q) = 0(A - 2M W ) — In 
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(2.16) 
(2.17) 



2q q 2 + 2AM n 

The N 3 LO corrections to the 2PE potential V$ have been recently calculated by Kaiser [18]. They 
arise from the one-loop "bubble" diagrams with both dimension two ttttNN vertices of the ci j ... j 4-type 
and from the diagrams which contain the third order pion-nucleon amplitude and lead to one-loop 
and two-loop graphs. We begin with the first group of corrections, for which one finds: 



167T 2 Fi 
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(218) 



No closed expressions can be given for some of the corrections from the second group.* 7 It appears 
to be convenient to give the contributions to the potential using the (subtracted) spectral function 
representation: 
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For the spectral functions pj(/x) (Vii^)) one finds: [18]: 
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(2.19) 



* 7 Entcm and Machleidt were able to calculate most of the integrals in eqs. (2.20) analytically (for A = oo) and to 
express the corresponding contributions to the potential in terms of the loop functions L°° (q) and A°° (q) [26] . 
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We use the scale-independent LECs d±, cfo, ds, ds, di4 and c/15 defined in [23]. In addition, one has 
to take into account the leading relativistic 1/m-corrections to the 2PE potential: 
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(2.23) 



Notice that these relativistic corrections differ from the ones given in ref. [14]. In fact, the specific 
form of the terms in eq. (2.23) depends on the form of the Schrodinger (or Lippmann-Schwinger) 
equation, see [25] for more details. The relativistic corrections given in eq. (2.23) are consistent with 
the relativistic Schrodinger equation (A. 4) and with the l/m 2 -corrections to the 1PE potential in 
eq. (2.13). 

Three-pion exchange starts to contribute at N 3 LO. The corresponding expressions for the spectral 
functions and the potential (obtained using dimensional regularization) have been given by Kaiser 
in [19, 20, 27]. It has been pointed out in these references that the 3PE potential is much weaker 
than the N 3 LO 2PE contributions at physically interesting distances r > 1 fm. Having the explicit 
expressions for the 3PE spectral functions, it is easy to calculate the potential in the SFR scheme. 
It is obvious even without performing the explicit calculations that the finite-range part of the 3PE 
potential in the SFR scheme is strongly suppressed at intermediate and short distances compared 
to the result obtained using DR. This is because the short range components which dominate the 
3PE spectrum are explicitly excluded in this approach. To illustrate that let us consider the isoscalar 
spin-spin contribution proportional to g\, which has been found in [20] to provide the strongest 3PE 
potential for 0.6 f m < r < 1.4 fm. The corresponding DR spectral functions Ps R 3lT (p) and p^^ip) 
are given by [20]: 



dr (lA g\(p-3M n f 

PS,3n{P) - 



35vr(32F3) 2 



DR , v _ 0a(M-3M w ) 
Pt,3tt{P) — 



M 3 /If 4 M 5 
2M 2 - 12/xM^ - 2 M 2 + 15^ + 2^f + 3^f 

\x p z p 6 

o , , „ Mi Ml Ml 



. (2.24) 



The finite-range part of the spin-spin potential Vs(r) can be obtained from the corresponding spectral 
functions via 

1 f 00 

V s( r ) = -7TT \ dwe-^ {p 2 PT (p) - 3 P s(p)) ■ (2.25) 
Using the spectral functions in eq. (2.24) one then gets for Vp^ n (r): 

4 — 3M r 

VsXir) = 2(8 ^ 2)3 e / (1 + ^) 2 (2 + M^rf . (2.26) 

Introducing the cut-off A in the spectral functions via p| F |^(/w) = Ps^ n (p) 0(A — p) and pj^^ n {p) = 
Pt%(p) #(A - n) we obtain for the SFR potential V§ F ^(r): 



120 + 120y + 60y 2 + 20y d + 5y 4 + y b (2.27) 



30(32vrF 2 ) 3 r 7 

- 25M 2 r 2 (6 + 6y + 3y 2 + y 3 ) + 45M 3 r 3 (2 + 2y + y 2 ) + 30M> 4 (1 + y) - 63M> 5 

where we have introduced the abbreviation y = Ar. In Fig. 1 we plot V^ v (r) and V§ F ^(r) for r from 
0.5 to 2 fm. The potential calculated using the spectral function regularization is much smaller in 
magnitude compared to the one obtained using dimensional regularization. Clearly, such a suppression 
does not take place at very large r, where Vg F ^(r) approaches Vp^ n (r). At such distances, however, 
the 3PE potential becomes negligibly small compared to the 1PE and 2PE contributions simply due 
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to its shorter range. As a consequence, the 3PE potential can be neglected everywhere except the 
region of very small r, where it anyhow becomes unreliable. This is further exemplified in Fig. 2, 
where we show the ratio of the N 3 LO isoscalar spin-spin contributions of 3PE and 2PE using both 
regularization schemes for a wide range of r. It turns out that the 3PE contribution reaches for r > 0.5 
fm at most 2% — 8% of the corresponding N 3 LO 2PE contribution depending on the choice of the 
spectral function cut-off. We therefore neglect all 3PE contributions in the present analysis. 

Although we have regularized the 2PE contributions by cutting off the large-mass components in 
the spectrum (or, equivalently, by explicitly shifting the corresponding short-distance components to 
contact terms), the resulting potential still behaves incorrectly at large momenta (or equivalently at 
short distances). The effective potential is valid for small values of the momentum transfer q and 
becomes meaningless for momenta q > A x . Moreover, since the potential V grows with increasing 
momenta q, the scattering equation is ultraviolet divergent and needs to be regularized. Following the 
standard procedure, see e.g. [16], we introduce an additional cut-off in the LS equation by multiplying 
the potential V(p, p') with a regulator function / A , 

V(p, p') - f\p) V(p, p>) / V) • (2-28) 

In what follows, we use the exponential regulator function 

f A (p) = exp[-p 6 /A 6 ] . (2.29) 

We will specify the values of the cut-offs below. 

It should be understood that our treatment of the effective potential is based on the heavy baryon 
formalism. As demonstrated in [28, 29], heavy baryon expansion becomes formally invalid for cer- 
tain two-pion exchange contributions at very large distances. This problem with the heavy baryon 
formalism has been first observed in the single-nucleon sector and can be dealt with using e.g. the 
Lorentz invariant scheme proposed by Becher and Leutwyler [30]. It is clear, however, that the NN 
interaction due to two-pion exchange becomes very weak at large distances, so that the problem with 
the formal inconsistency of the heavy baryon approach is expected to have little relevance for practical 
applications. 

Last but not least, we would like to comment on some key features of the SFR scheme adopted in the 
present work. First of all, it is crucial to understand that this approach does not affect the "chiral 
features" of the NN potential. We remind the reader that the resulting effective NN potential consists 
of the long- and short-range pieces. Spontaneously broken approximate chiral symmetry of QCD leads 
to highly non-trivial constraints for the the long-range part of the potential, which is given by the 
terms nonpolynomial in momenta. The short-range part of the potential given by a series of the most 
general contact interactions with increasing power of momenta is not affected by chiral symmetry with 
the exception of the quark-mass dependence of the corresponding LECs, which is not relevant for the 
present analysis. In other words, only the long-distance asymptotics of the potential is constrained by 
chiral symmetry. The SFR scheme does, per constraction, not affect the long-distance asymptotics of 
the potential and leads to the same result as obtained using DR, see [17] for more details. The only 
difference to the DR result is given by a series of the short-range interactions. It is, therefore, obvious, 
that the the SFR method does not affect the constraints of the chiral symmetry implemented in the NN 
potential. Further, we point out that the equivalence of the SFR and the finite cut-off regularization 
has only been established at a one-loop level and does not hold true for both loop integrals of the 
N 3 LO 2PE contribution. The prominent feature of the applied regularization scheme is given by the 
fact, that it only affects the two-nucleon interaction. One can, therefore, directly adopt the values 
for various LECs resulting from the single-nucleon sector analyses, where dimensional regularization 
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has been used. On the contrary, if a finite momentum cut-off regularization would be applied to 
both loop integrals entering the N 3 LO 2PE contribution, one would need to re-extract the values of 
the corresponding LECs from pion-nucleon scattering and the process ttN — ► 7T7riV using the same 
regularization scheme. 



2.4 Isospin breaking effects 

Isospin-breaking nuclear forces have been extensively studied within effective field theory approaches, 
see e.g. [31, 32, 33, 34, 7, 35], as well as using more phenomenological methods, see e.g. [36, 37] for 
some recent references. In the Standard Model, isospin-violating effects have their origin in both 
strong (i.e. due to the different masses of the up and down quarks) and electromagnetic interactions 
(due to different charges of the up and down quarks). The electromagnetic effects can be separated 
into the ones due to soft and hard photons. While effects of hard photons are incorporated in effective 
field theory by inclusion of electromagnetic short distance operators in the effective Lagrangian, soft 
photons have to be taken into account explicitly. 

In the present analysis we are rather limited in the treatment of isospin-violating interaction, which 
have to be included precisely in the way it is done by the Nijmegen group [38]. This is due to the fact 
that we are using the Nijmegen phase shifts instead of the real data as an input to fit the unknown 
LECs. Let us explain this point in more detail. With the only exception of the 1 Sq partial wave, the np 
isovector phase shifts in the Nijmegen PWA are not obtained independently from np data, but rather 
extracted from the proton-proton (pp) phase shifts using the assumption that the differences in the 
phase shifts result entirely due to isospin-breaking effects associated with m p / m n and M n ± / M n o 
in the 1PE potential as well as due to electromagnetic interactions. In order to be consistent with the 
Nijmegen phase shift analysis, we therefore have to neglect various isospin-breaking corrections and 
adopt the same isospin-breaking and electromagnetic interactions as in [38]. Nevertheless, we have 
decided to overview the dominant isospin-breaking contributions and to remind the reader on their 
relative size following mainly the lines of reference [7] but extending the consideration to higher orders. 
For a detailed review of charge-symmetry breaking in the nucleon-nucleon interaction the reader is 
referred to [35]. 

Consider first isospin breaking in the strong interaction. The QCD quark mass term can be expressed 
as 

C^ = -^q(m Q + m d )(l-er 3 )q , (2.30) 

where 

fnd — mu 1 

e = — ~ - . 2.31 

m d + m u 3 

The above numerical estimation is based on the light quark mass values utilizing a modified MS 
subtraction scheme at a renormalization scale of 1 GeV. The isoscalar term in eq. (2.30) breaks chiral 
but preserves isospin symmetry. It leads to the nonvanishing pion mass, M% = {mu + m^B / 0, where 
B is a low-energy constant that describes the strength of the bilinear light quark condensates. All 
chiral symmetry breaking interactions in the effective Lagrangian are proportional to positive powers 
of M%. The isovector term (oc T3) in eq. (2.30) breaks isospin symmetry and generates a series of 
isospin-breaking effective interactions oc (eM^) n with n > 1. It therefore appears to be natural to 
count strong isospin violation in terms of eM%. However, we note already here that isospin-breaking 
effects are in general much smaller than indicated by the numerical value of e, because the relevant 
scale for the isospin-conserving contributions is A x rather than m u + ma- 

Electromagnetic terms in the effective Lagrangian can be generated using the method of external 
sources, see e.g. [39, 40, 41] for more details. All such terms are proportional to the nucleon charge 



12 



matrix Q c ^ = e(l + t^)/2, where e denotes the electric charge.* 8 More precisely, the vertices which 

contain (do not contain) the photon fields are proportional to Q™ h (Q^), where n = 1,2, Since 

we are interested here in nucleon-nucleon scattering in the absence of external fields, so that no 
photon can leave a Feynman diagram, it is convenient to introduce the small parameter e 2 ~ 1/10 for 
isospin-violating effects caused by the electromagnetic interactions. 

Due to its perturbative nature induced by the small parameters eM 2 and e 2 , we treat the strong and 
electromagnetic isospin violation in addition to the power counting of the isospin symmetric potential 
mentioned in section 2.1. Although not necessary, in practical applications it often appears to be 
more convenient to have a single expansion parameter. Thus, one has to relate the quantities e, e to 
the generic low-momentum scale Q related to external three-momenta of nucleons and the pion mass 
(p ~ p' ~ M w ~ Q) and introduced before. Here and below, we will make use of the following simple 
counting rules* 9 : 

Q 

e ~ e ~ -p . (2.32) 

The counting of the electric charge is consistent with the one commonly used in the pion and pion- 
nucleon sectors, see e.g. [39, 42, 43, 44] (it differs, however, from what is commonly used in the 
description of extremely non-relativistic hadronic bound state, see e.g. [45]). In addition to the above 
mentioned counting rules, we need to deal with the extra l/(47r) 2 -factors, which typically arise when 
calculating loop integrals. For pion loops, such factors are naturally incorporated in the chiral power 
counting through the relation A x ~ 4irF n . For photon loops we will further assume, that 

e 2 O 4 

W (2.33) 



(4vr)2 A4 ' 

which simply means that the factors l/(47r) 2 provide two additional powers of the small parameter. 
In the following, we will denote the order of various isospin-violating interactions by "L0", "NL0', 
... in order to distinguish the above mentioned phenomenological extension of the counting rules from 
the usual chiral power counting in the isospin-conserved case. Certainly, one has always the option 
to discard this generalization of the chiral counting rules and to perform separate expansions in e, 
e and Q/A x . Notice further that eqs. (2.32) and (2.33) suggest a different counting of the strong 
isospin-breaking effects compared to ref. [7]. In that work strong and electromagnetic effects have 
been classified using a separate expansion without introducing a unified expansion scheme. The L0, 
NL0, . . . contributions in the present work should therefore not be confused with the corresponding 
terms in [7]. Last but not least, the above counting scheme is similar to the one adopted e.g. in 
[31, 33], where effects ~ a/ir were also considered as being one order suppressed compared to the ones 
~ ^M 2 /A 2 . 

Let us now apply the power counting rules to estimate isospin-violating corrections to hadronic masses, 
see also [35] for a similar estimation. We begin with the pion mass. It is well known that the pion 
mass does not receive contributions linear in the quark mass difference and the strong contribution to 
the pion mass starts at the second order in rrid — m u . Consequently, the leading strong term can be 
estimated as (AM 2 ) str = (M 2 ± - M 2 ) str oc (eM%) 2 A~ 2 , which is of the order v = 6 to be compared 
to v = 2 for the isospin-symmetric term oc M 2 . One thus expects the strong contribution to the 
pion mass difference (AM 7r ) str = (M n ± — M^)) str to have the size ~ 0.1 . . .0.3 MeV depending on 

* s Or equivalently, one can use the quark charge matrix e (1/3 + ts)/2. 

* 9 This suggests a slightly different counting of the strong isospin-breaking effects as compared to [7]. Most of the 
conclusions of [7] remain, however, unchanged. The important difference is that the leading isospin-violating short-range 
interaction is now proportional to the quark mass difference, while electromagnetic contact terms are shifted to higher 
orders. 
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whether one substitutes in the numerical estimation M p or 4irF n for A x . Clearly, we cannot predict 
whether the shift is positive or negative. The leading electromagnetic contribution to AM 2 is of the 
order (AM 2 ) em ~ e 2 (47r) _2 A 2 , which is a v = 4-effect according to our counting rules. Numerically, 
one estimates the size of the pion mass difference to be (AM n ) 

cm ~ 1 ... 3 MeV. We see that both 
power counting arguments and numerical estimations suggest that the pion mass difference is mainly 
of electromagnetic origin. Furthermore, the estimated size of the electromagnetic shift agrees well with 
the observed value M^i-M^o = 4.6 MeV. All these statements can also be backed by hard calculations, 
for a classical review see [46]. For the nucleon mass difference, the strong contribution is linear in 
the quark mass differences and can be estimated as: (Am) str = (m n — m p ) st i ~ (eM^A" 1 ~ 6 ... 9 
MeV.* 10 According to the counting rules, this is the v = 3-effect. Electromagnetic shift appears at 



v = 4 and is expected to be of the order (Am) 



: (4tt)" 2 A x ~ 0.5 . . . 0.7 MeV. In reality, the effects 



are of opposite sign and the difference between them is less pronounced. One observes (jn n 



m 



pjstr — 



2.1 MeV and (m n — m p ) em 
again [46] for more details). 



~ —0.8 MeV, leading to the physical value of m n 



m. 



1.3 MeV (see 



We are now in the position to discuss various isospin-breaking contributions to the two-nucleon force. 
As explained in [7], the leading-order (i.e. L0) isospin-breaking interactions are due to the pion 
mass difference in the 1PE potential and the static Coulomb interaction. The latter is clearly of the 



order* 11 



2 F%, while the former is 



AM 2 
M 2 



^— A 2 

(4vr) 2 x 



Ml 



O 



91 



(2.34) 



where we used of the counting rules (2.32) and (2.33) together with M w ~ Q. Thus, the L0 isospin- 
breaking force is of the order v = 2. Consider now NL0 corrections to this result, which appear 
at v = 3. The pion-nucleon coupling constant receives strong isospin- violating contributions of the 
order eM 2 /A 2 . The corresponding LECs in the pion-nucleon Lagrangian are denoted by dn, d\s 
and dig in the notation of Ref. [23] . This charge dependence of the pion-nucleon coupling constant 
leads to isospin-violating 1PE of the order v = 3. In addition, one has to take into account strong 
isospin-breaking contact interaction of the kind 



eM 2 (N^t 3 N) (N^N) 



(2.35) 



which leads to charge symmetry breaking. We can check the accuracy of our estimation numerically 
using the values for the LECs found in [7]. According to eq. (2.35), we expect the ratio of the isospin- 
breaking terms to isospin-conserving ones to be typically of the size: eM 2 /A 2 ~ 0.5% . . . 1.1%, where 
the uncertainty results again from using two different estimations for A x . Picking up the numbers 
from Table 2 in [7] we find for this ratio the values 0.8%, 0.8% and 3.1% for three different values of 
the (sharp) cut-off A in the Lippmann-Schwinger equation: A = 300, A = 400 and A = 500 MeV. 
Thus, our numerical estimation is consistent with the results of [7]. The NNL0 corrections are of the 
order v = 4 and arise from various sources. First, one has to take into account isospin-breaking in 
the 2PE potential due to electromagnetic corrections to the pion-nucleon coupling, see e.g. [33]. The 
correction due to the pion mass difference in the leading 2PE potential can be estimated as: 



AM 2 Q 2 _ 
Ml A 2 



Q 



ii 



A 4 

x 



(2.36) 



* 10 This too large value reflects our earlier statement about the use of the parameter e to estimate isospin-breaking 
corrections. 

* n The factor F% results from the common normalization of the isospin-symmetric part of the two-nucleon potential 
adopted in this work. This factor can be understood e.g. from looking at the 1PE potential in eq. (2.9). 
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Another isospin-violating two-pion exchange interaction at this order is generated by the triangle and 
football diagrams with one insertion of the isospin-breaking irirNN vertex with the LEC C5 (using 
again the notation of [23] ) . This vertex is proportional to eM 2 and is thus formally of the lower order 
than the electromagnetic ttttNN vertices oc e 2 /(4-7r) 2 . As we have seen on an example of the nucleon 
mass difference, in practice, both effects might be of a comparable size. For a recent work on this kind 
of isospin-breaking forces see [37]. The remaining contributions are given by the static ^-exchange 
of the order ~ e 2 /(47r) 2 and by two independent contact interactions 



e 2 ,„ + + _ . e 2 



, , 2 (iVV 3 iV) (iVVsiV) and (AV3A) (A f N) (2.37) 

which lead to both charge independence and charge symmetry breaking. Notice that, in practice, the 
effect of the second interaction cannot be disentangled from the effect of term in eq. (2.35). 

Let us now estimate the size of isospin violation in the NN scattering due to the nucleon mass difference. 
We first note that the first relativistic corrections to the isospin-symmetric part of the two-nucleon 
force appear at at N 3 LO [y = 4) and are given by l/m 2 -corrections to the 1PE and 1/m-corrections 
to the leading 2PE potential. Consequently, the size of the corresponding isospin-violating terms can 
be estimated as 



Am Q 4 _ eM 2 Q 4 = Q 



m A£ mA x A 4 



(2.38) 



Such terms therefore contribute only at the order v = 8. In addition to the above mentioned correc- 
tions, one has to account for the fact that the neutron-proton mass difference leads to energy shifts of 
virtual states when calculating two-pion exchange diagrams. This can also easily been understood in 
the language of the heavy baryon formalism: factoring out the exponential factor exp(im p v • x) from 
the proton and neutron fields, where v and x denote the proton velocity and position, the neutron 
propagator receives db shift in the denominator oc {rn n — m p ) after integrating out the small field compo- 
nents. It is then easy to see that the isospin-violating 2PE is suppressed against its isospin-conserving 
part by a factor: 



Am eM 2 



r 2 l 

(2.39) 



Q A X Q ~ LA 2 

Therefore, neutron-proton mass difference in 2PE starts to contribute at v = 4. This sort of charge 
symmetry breaking corrections has been studied recently in [36, 35]. Notice further that, as pointed 
out in [35], certain loop integrals in the 2PE contributions give only one power of (Air) instead of 
expected two powers and are, therefore, enhanced. We will not take this enhancement into account 
in the present work. Apart from the above mentioned corrections to the nucleon-nucleon force, the 
neutron-proton mass difference has to be taken into account in kinematical relations as discussed in 
appendix A, as well as in the expression for the kinetic energy of the nucleons. Let us consider this 
last effect. Its contribution to the scattering amplitude can be estimated by looking at the Lippmann- 
Schwinger equation 

T = V + VG T , (2.40) 

where Go refers to the free propagator of two nucleons. Both terms on the right-hand-side of the 
above equation are of the same order v = 0. Taking into account the nucleon mass difference in Go 
leads therefore to a correction to the T-matrix of the order 



Am eM 2 

~ — » — = ^ 

m mAv 



At 



(2.41) 

and thus contribute at v = 4. Notice that this estimation is valid for both relativistic and nonrela- 
tivistic expressions for the two-nucleon propagator Gq. 
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All other isospin-violating corrections are suppressed by further powers of the small parameter. We 
would like, however, to point out an important limitation of our estimation due to the fact that we 
do not explicitly account for the long-range nature of electromagnetic forces. Consider, for example, 
the leading one-pion and one-photon exchange forces. For simplicity, we will restrict ourselves to the 
1 Sq proton-proton channel, where the 1PE potential takes the form: 



The static Coulomb interaction, 



e 2 



Vcoulomb(<Z) = ^ > ( 2 - 43 ) 

is suppressed compared to vf® (q) by two powers of the small parameter Q/A x according to the power 
counting. Such an estimation works fairly well for momenta q of the order q ~ M n , for which we 
get V^(q = M n ) ~ 23 GeV~ 2 and Vc uiomb((7 = M w ) ~ 5 GeV~ 2 . The power counting, however, 
breaks down for small momenta q <C due to the long-range nature of the Coulomb interaction. 
For example, for q = M n /4 one gets: V$(q = M w /A) ~ 3 GeV~ 2 while V Co uiomb(<7 = M n /A) ~ 82 
GeV -2 . Consequently, the Coulomb interaction provides the dominant contribution to the potential 
for small momenta and requires a nonperturbative treatment at low energy. A possible way out of 
the above mentioned inconsistency would be to develop separate and systematic power counting for 
momenta much smaller than the pion mass. This is, however, beyond the scope of the present work. 
Notice that a similar idea with two different power counting regimes has been applied recently to the 
nucleon Compton scattering in order to extend the region of applicability of the effective field theory 
in the A-region [47]. In the present analysis, we will simply take into account higher-order corrections 
to the long-range electromagnetic interactions when determining the values of the LECs in order to 
correct for the low-momentum behavior of the NN potential. The first long-range corrections beyond 
the ones considered above result from two-photon exchange, whose size can be estimated as 



p 2 p 2 



Q 2 (4vr) 2 77 



91 



(2.44) 



It thus formally appears at the order v = 6. In addition, at the same order v = 6 one has to take into 
account relativistic 1/m corrections to the static one-photon exchange, which provide a contribution 
of the following size: 



i_ F 2 9- = JL F 2 = o 

Q2 r n m 2 m 2 r n 



Q 



A 6 

x 



6" 

(2.45) 



The relative sizes of various isospin-breaking contributions discussed above are summarized in Table 1 . 
In what follows, we will give explicit expressions for the above mentioned interactions. 

2.4.1 Finite range isospin breaking forces 

Let us now give the explicit expressions for the finite-range isospin-violating interactions up to NL0. 
The dominant v = 2 contribution (L0) due to M w ± / M n o can be taken into account by replacing the 
isospin-conserving expression V\- K (q) in eq. (2.13) by: 

9a\ 2 -. -*oi-qo 2 -q(, p 2 + p' 2 



v f\ \r < \ ( 9A \ -,B\-qBi-q ( 



2m 2 



16 



l7r, np, T- 



-o(q) 



9A 
2F„ 



ri ■ r 2 



< ?i ■ q <?2 ■ q 

q 2 + M 2 



1 



p 2 + p' 2 
2m 2 



(2.46) 



where T denotes the total isospin, M^± and M n o are the masses of the charged and neutral pions, 
respectively, and 

Mir = + = 138.03 MeV . (2.47) 



The 1PE potential gets further charge independence and charge symmetry breaking contributions at 
NL0 and NNL0 due to isospin violating pion-nucleon couplings. The final expression for the 1PE 

potential is then of the kind: r\iV^ , where rji are the channel-dependent constants: rj nn / i] np / r] pp . 
Unfortunately, the actual size of isospin-violating corrections to the pion-nucleon coupling is not well 
determined at presence [33]. 

The pion mass difference in the 2PE can be incorporated as outlined in ref.[48]. It is most convenient 
to consider the isoscalar and isovector 2PE piece separately, 



V 27 r = Vl + Vln-T 2 . 

The isoscalar part V 2 \ can be expressed as [48] 

Vi = 2 -Vl(M^,M^ + ) + \vl{M %0 , M n o) = M n ) + O 



(2.48) 



(2.49) 



where the arguments of V 27r denote the masses of exchanged pions. For the isovector 2PE V 2 \, one 
has the general structure 



V, 



2tt 



t? r| V 2 \{M„ + ,M„ + ) + (n • r 2 - rf r|) V 2 1 7T (M 7T + , M w o) 
V 2 \{M„ + ,M„ + ) 



for pp and nn , 



2V 2 \ (M w+ , M w0 ) - V 2 \ (M w+ , M w+ ) ~ V 2 \ (M w0 , M T o ) for np, T = 1 . 



(2.50) 



The result in the last line of the above equation is valid modulo ((M n + — M 7r o)/M 7r +) 2 -corrections. 
For the T = case the 2PE potential reads V 27V (M 7V , M n ). 

The 7T7 exchange diagrams have been calculated in ref.[32] and we give below the results obtained 
in that paper omitting all computational details. Due to isospin, only charged pion exchange can 
contribute to the irj potential V ni and thus it only affects the np system. The potential has the form 



K-7O7) 
K 7 (/3) 



AF 2 M 2 ( f i ' ^2 - Ti 72)01 • qd 2 -qV^{P) , 
2/? 4 (l + (3 2 ) { P ' 2/3 2 1 + (5 2 



(2.51) 



Here, = \q\/M n + and 7 is a regularization scheme dependent constant. The analytical form of V nJ 
is similar to the one of the 1PE potential, but it differs in strength by the factor a /-k ~ 1/400. 

Finally, the expressions for the remaining isospin-violating 2PE contributions at NNL0 have been 
given in [35]. 
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2.4.2 Long range (soft) isospin breaking forces 

We now discuss long-range isospin-breaking interactions which are often referred to in the literature 
as "electromagnetic forces". The static Coulomb force in eq. (2.43) does certainly not completely 
represent the electromagnetic interaction between two nucleons but only its leading contribution. 
The first long-range corrections to the static Coulomb force are either suppressed by m~ 2 (relativistic 
corrections to the static one-photon exchange) or by an additional power of the fine-structure constant 
a (two-photon exchange). Although all these effects are formally of higher order, we nevertheless 
prefer to take them into account explicitly for the following reasons. First of all, the effects of these 
interactions are magnified at low energy due to their long-range nature. Further, as explained above, 
in our analysis we have to take into account isospin-breaking effects in the same way as it is done in 
[38] . The electromagnetic interaction for the pp and np case is given by 



Vem(pp) = Vci + Vc2 + Vyp + Vmm(pp) , 
V E M(np) = V M M(np) , 
V E M(nn) = V M M(nn) , 



(2.52) 



where Vci and Vq2 are usually referred to as "improved Coulomb potential" . They take into account 
the relativistic l/m 2 -corrections to the static Coulomb potential and include contributions of the 
two-photon-exchange diagrams [49]. The explicit coordinate-space expressions read: 



Vci = 
V C 2 = 



a 
r 



2m 2 p 



(A + ^ + ^(A + fc 2 ) 



aa 



(2.53) 



where A denotes the Laplacian. The energy-dependent constant a' is given by 



a 



a- 



m 2 p + 2k 2 



m p J m 2 + k 2 



(2.54) 



Here k is the c.m.s. scattering momentum. The term Vq2 is chosen in such a way that it leads to an 
exact cancellation between the proper two-photon and the iterated one-photon exchange, see [49] for 
more details. The approximation made in the second line of eq. (2.53) is based upon using Coulomb 
distorted-wave Born approximation (CDWBA), see [50] for more details. The modified Coulomb 
potential Vqi in eq. (2.53) can be treated in momentum space in the same way as the usual static 
Coulomb potential as described in appendix C.* 12 The magnetic moment interaction Vmm in eq. (2.52) 
is given by [51]: 



Vmm(pp) 



Vmm (ran) 



a 



4 m 2 r 3 



V M M.(np) = - 



CtKr, 



2m„r 3 
4m 2 r 3 



vlSn + (6 + 8k p )L ■ S 



^ S 12 + -(L.S + L.A) 



2m 



m 



S 



12 ■ 



(2.55) 



where fi p = 2.793 and fi n = —1.913 are the proton and neutron magnetic moments and n p = fi p — 1, 
K-n = fJ>n their anomalous magnetic moments. Further, L is the orbital angular momentum, A = 

#12 Clearly, one has to use the appropriately adjusted regular and irregular Coulomb functions Fi(r) and Gi(r). 
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(c?i — 02)12 and S^i = (<?i ■ r)((?2 ■ r)/r 2 — o\ ■ (J 2 /3. The corresponding expressions in momentum 
space can be found e.g. in [52]. Finally, the vacuum polarization potential Vyp derived by Ueling [53], 
see also [54], reads: 

where m e is the electron mass. Clearly, the vacuum polarization potential is not of an infinitely long 
range. Its range is governed by the electron mass, which is still tiny compared to the relevant mass 
scales in the nucleon-nucleon problem. This is similar to the treatment of vacuum polarization in 
EFT approaches for hadronic bound states, see e.g. [55]. 

It is important to realize that the expressions (2.52) refer to point-like nucleons and only define 
the long-distance asymptotics of the corresponding electromagnetic interaction. The short-distance 
structure is more complicated and not shown explicitly. In particular, we do not include zero-range 
(for point-like nucleons) terms as well as electromagnetic form factors which can, in principle, be 
calculated consistently in EFT. Such short-range terms with the nucleon form factors of a dipole 
form are, for example, included in the Argonne V18 potential. Last but not least, we note that 
the above consideration of the electromagnetic effects is based on the "nonrelativistic" Schrodinger 
equation (3.12), which will be defined in the next section. To close this section let us point out some 
well known practical complications which arise due to the presence of the long-range electromagnetic 
forces. 

• Asymptotic states are affected by electromagnetic interactions. The S-matrix has to be formu- 
lated in terms of asymptotic Coulomb states. 

• The formally suppressed (as compared to the strong nuclear force) electromagnetic interactions 
are enhanced at low energy. The Coulomb interaction requires a nonperturbative treatment. 
Even the effects due to magnetic moment interaction might be large for certain observables 
under specific kinematical conditions. For example, in the np system, it gives rise to a forward- 
angle dip structure for the analyzing power. 

• The expansion of the scattering amplitude in partial waves converges very slowly in the presence 
of magnetic moment interactions. 

3 Scattering equations 

We start with the relativistic Schrodinger equation (A. 4) and assume the potential to be of a fi- 
nite range. The treatment of the nucleon-nucleon scattering problem in the presence of the long- 
range Coulomb interaction is relegated to appendix C. The scattering states are described by the 
Lippmann-Schwinger equation corresponding to the Schrodinger equation (A. 4). The LS equation 
(for the T-matrix) projected onto states with orbital angular momentum /, total spin s and total 
angular momentum j is 

v , Jo \ Z7r ) 2\Jp A + m z — 2 VP +m z + ir] 

(3-1) 

with rj — ► + . In the uncoupled case, I is conserved. The partial wave projected potential Vjfj(p',p) 
can be obtained using the formulae collected in appendix B. The relation between the S— and on- 
the-energy shell T-matrix is given by 

W = ^ - ^Py/F+^Tg(p) . (3.2) 
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The phase shifts in the uncoupled cases can be obtained from the 5-matrix via 

S% = e 2i ^\ S}> = <?»¥, (3.3) 



where we have used the notation . Throughout, we use the so-called Stapp parametrization [56] of 
the S-matrix in the coupled channels (j > 0): 

5=| S i- 1 J- 1 S i~ U • 1 | 



i(2t 

i sin (2e) exp (^j^ + ^j+i) cos (2e) exp (2i<5j^ 



cos (2e) exp (2i(5^'£_ 1 ) i sin (2e) exp (i^ij + ^j+i) 



(3.4) 



For the discussion of the effective range expansion for the 3 Si partial wave we will use the different 
parametrization of the 5-matrix, namely the one due to Blatt and Biedenharn [57]. The connection 
between these two sets of parameter is given by the following equations: 

<5j_i + Sj + i = <5j_i + 5j + i , 
tan(2e) 

sm^-^i) = 
•fx x ^ sin(2e) 

sm^-^i) = -^i), (3.5) 

where 5 and e denote the quantities in the Blatt-Biedenharn parametrization and we have omitted 
the superscripts for 5's. 

To close this section we would like to remind the reader that the Schrodinger and Lippmann-Schwinger 
equations (A. 4) and (3.1) may be cast into a nonrelativistic form. One way to do that is using the 
Kamada-Glockle transformation [58], which relates the relativistic and nonrelativistic c.m.s. momenta 
p and q via: 

T kin = 2vV + m 2 - 2m = ^- . (3.6) 

m 

The potential V to be used in the nonrelativistic Schrodinger equation 



r g 2 

1 



m 

is defined in the partial-wave projected representation as 



4> = E4>, (3.7) 



-.1 2 



4m 2 ' 
(3.8) 

where V^{p, p') is the potential entering the relativistic Schrodinger equation (A. 4). The wave- 
function 4> is related to \I> in eq. (A. 4) via 



2 r A 



«"=vi i+ 2&)v i+ 4&*(y< 2+ 4&j- < 3 - 9 » 

The S-matrix is defined via 

^(Q)=Sw--^Qmf^( q ), (3.10) 
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where the T-matrix TfJ satisfies the usual nonrelativistic Lippmann-Schwinger equation 

f*jM) = i*W)+x; I %^y^') ^_™ 2+i f^\<i) . (3.ii) 

It can be demonstrated [58] that the S-matrix S^, equals for any given energy the S-matrix S^, defined 
in eq. (3.2), that is <§^(g) = S^,{p). Another commonly used way to cast the relativistic Schrodinger 
equation (A. 4) into a nonrelativistic-like form is based upon the algebraic manipulations with this 
equation, see [25]. More precisely, adding 2m to both sides in eq. (A. 4) with subsequent squaring 
them, subtracting 4m 2 and dividing both sides by 4m leads to 



— + V 



k 2 

= — (3.12) 
m 



m 

where the momentum k is related to the energy E in eq. (A. 4) via 

E = 2^k 2 + m 2 -2m, (3.13) 
and the potential operator V is given by 

or, in the partial-wave projected basis, by 

i>sj, f\ ( vV + m2 , vV 2 + m2 \ rrsj, f\ . 1 sr> f°° dp"p" 2 sj „ sj „ , 
v u ?(p,p ) = [-^^ + 2m j V(p.p ) + ^ E y o V(p.p WAp >p ) ■ 

(3.15) 

The curly bracket in eq. (3.14) denote an anticommutator. Notice that contrary to the previously 
described approach, the "nonrelativistic" Schrodinger equation (3.12) still requires relativistic kine- 
matics in relating the energy and momentum, see eq. (3.13). The S- and T-matrices Sfj 1 , and T^, 
are defined via eqs.(3.10) and (3.11), respectively (with S, T, q, q' being replaced by S, T, p, p'). At 
any given momentum p one has S^, (p) = S 1 ^ (p) . We have also checked numerically that both equa- 
tions (3.7) and (3.12) lead to identical results. It should be understood that both ways to cast the 
relativistic Schrodinger equation (A. 4) into a nonrelativistic form discussed in this section are limited 
to the two-nucleon problem. To the best of our knowledge, no extension to different systems has yet 
been offered. Consequently, three- and more-nucleon observables calculated using a nonrelativistic 
approach with the NN potential V or V will lead to different results. One should therefore use the 
relativistic Schrodinger equation (or Faddeev-Yakubovsky equations) with the potential V in such 
cases. The same applies for processes with external probes. 

4 Bound state 

We now turn to bound state (i.e. deuteron) properties. The deuteron binding energy E<\ and wave 
function tyf (p) can be obtained from the homogeneous part of eq.(3.1): 

*?(p) = ; 1 T r tS^t VitM) *?'(p') , (4- 1 ) 

' E d - (2y/p 2 + m 2 - 2m) y Jo ( 27r ) 3 ' 
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with s = j = 1 and I = I' = 0, 2, or, alternatively, from the nonrelativistic-like equation 



*?(p) 



£ d + £^/(4m) - p 2 /m ^ J (2vr 



E 



(4.2) 



where t^y,(p,p') is related to V^ip^p') via eq. (3.15). Here we have used the relation (3.13) between 
the binding energy and momentum. In addition, one can also use the nonrelativistic Schrodinger 
approach as described in the previous section, which leads to: 



4(p) 



E 



E d -p 2 /m^J (2vr) 



00 dp'p' 2 



(4.3) 



where V, S u(p,p') is defined via eq. (3.8). The wave functions 4>f(p) and *&f(p) are related via eq. (3.9). 

We will now regard the so-called static properties of the deuteron using the nonrelativistic-like equa- 
tion (4.2). The latter is fully equivalent to the relativistic equation (4.1) and leads to the same wave 
function, but has the advantage that one can apply the standard nonrelativistic formulae to study 
various deuteron properties. We denote by u(r) and w(r) the S- and D-wave components of the 
coordinate space wave function ^f(r) and by u(p) and w(p) the momentum space representations of 
u(r)/r and w(r)/r: 

2 f°° 2 f°° 

u (p) = — / u(r)jo(pr)rdr and w(p) = — / w(r)jo(pr)rdr . (4-4) 

K J0 ' 7T Jo 

The wave functions u and w are normalized according to: 

poo roo 

/ dpp 2 [u(p) 2 + w(p) 2 } = dr [u(rf + w(r) 2 } = 1 . (4.5) 
Jo Jo 

The probability Pd to find the nucleons inside of the deuteron in a D-state can be calculated via 

roo roo 

Pd = dpp 2 w(p) 2 = / drw(r) 2 . (4-6) 
Jo Jo 

Further, one can compute the deuteron quadrupole moment Qd and the matter root-mean-square 
(rms) radius y 7 (r 2 )^ through the following equations: 



1 f°° 

Qd = — / dr r 2 w(r) [y/8u(r) — w (r) 
20 Jo 



dp dp 



du(p) 
dp 



+ p 



( dw(p) V 
V dp J 



+ Qw{pY 



and 



,2\d 



roo 

I drr 2 [u(r) 2 + w(r) 2 ] 
Jo 



1/2 



The wave functions u(r) and w(r) behave at large r as: 



3 3 

u(r) ~ A s e~ yr , w(r) ~ A D e~' yr ( H h 



ryj, (7r) S 



(4.7) 



(4.8) 



(4.9) 
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where A$ and Ad are the asymptotic normalization factors of the S- and D-states, respectively, 



and 7 = \J\mE& + E 2 /A\. Instead of the quantities As and Ad, one often introduces the deuteron 
normalization and the asymptotic D/S ratio % according to: 

N 2 = A 2 +A 2 D , % = ^- (4.10) 

Not all of the above mentioned deuteron properties are observable and can be measured experimentally. 
The D-state probability is well known to be unobservable [59]. The deuteron electric quadrupole 
moment corresponds to the quadrupole form factor at \q\ = 0, where q denotes the momentum trans- 
fer. Clearly, the expression (4.7), which gives just the deuteron expectation value of the quadrupole 
operator Qij 

Qij^l^irj-Sijr 2 ), (4.11) 

is only an approximation to the experimentally measured value for the quadrupole moment, which 
i.e. does not take into account two-nucleon currents and relativistic corrections, see e.g. [60] for 
more details. A related discussion in the framework of EFT can be found in [61]. The situation is 
similar with the deuteron matter rms-radius y (r 2 )^, which is related to the experimentally measured 

deuteron charge rms-radius \J (r 2 )^ h via [62, 63, 64]: 

(r 2 ) d ch = (r% + (r 2 ) p ch + (rY ch , (4.12) 

where yj (r 2 ) p ch = 0.886(11) fm (taking the mean of the three recent values form Refs.[65, 66, 67] and 

adding the errors in quadrature) and (r 2 )™ h = —0.113(5) fm 2 [68] are the proton and neutron ms-radii, 
respectively, and the "point-nucleon" radius of the deuteron (r 2 )p t is given by 

(rY p t = ( r2 )i + (rY B - (4-13) 

Here (r 2 )^ subsumes the "nuclear" effects due to two-body currents as well as relativistic corrections. 
Notice that while the "point-nucleon" deuteron radius is measurable, the matter radius (r 2 )^ is clearly 
not an observable quantity. In particular, the separate contributions (r 2 )^ and (r 2 )^ change by a 
unitary transformation in the two-nucleon system, see [64] for more details. For one specific choice of 
such a transformation, the effects due to two-nucleon currents in (r 2 )^ are estimated to be of the order 
~ 0.016 fm 2 [64]. Contrary to the previously discussed deuteron quadrupole moment and rms-radius, 
the asymptotic quantities As and Ad (or, equivalently, and %) as well as the deuteron binding 
energy are observables related to the "pure" nucleon-nucleon system. In particular, the binding energy 
gives the position of the NN S-matrix pole, while the normalization is related to the residue of 
the pole in the following way, see e.g. [52, 69]: 

iV 2 = lim { -P-^l p mfll(p) . (4.14) 

Here fed = 27 and TJq(p) is the I = V = 0-component of the diagonalized T-matrix T 11 = UT XX U~ X , 
where 

U=( C ° S l Sin ^ V (4.15) 
y — sin e cos e J 

and e is the Blatt and Biedenharn mixing angle [57]. Alternatively, one can rewrite equation (4.14) in 
terms of the Blatt and Biedenharn eigenphase shift 5q(j>) as: 

iV 2 =hm 2l{p ~ kd) . (4.16) 
P -»fe d 1 _ ^tan [So(p)\ 
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Notice that the T-matrix becomes real and the phase shift Sq(p) imaginary at negative energies. 
Finally, we point out that the asymptotic D/S ratio rjd is given by the negative of the Blatt and 
Biedenharn mixing angle at the deuteron pole: [52]: 

?? d = -tane(/c d ) . (4-17) 



Up to now we have discussed the deuteron properties in the context of the nonrelativistic-like equation 
(4.2). As already pointed out before, one could alternatively use the nonrelativistic equation (4.3). 
Both schemes are completely equivalent for the two-nucleon system and lead to the same phase shifts 
and the deuteron binding energy. It is also clear from eq. (4.17) that the asymptotic D/S ratio rjd 
does not change when one uses eq. (4.3) instead of eq. (4.2). On the other hand, the normalization 
iVd or, equivalently, the asymptotic normalization As will change. Eq. (4.16) takes the form 

A> d 2 =lim 2 ^"! d) . (4.18) 
p^fe d i _ j tan [S (P)] 

if one uses the nonrelativistic Schrodinger equation (4.3). Here, = i^/\mE^\ is the nonrelativistic 

deuteron binding momentum and Sq is the S-wave eigenphase shift calculated using eq. (3.11). Since 
both schemes are phase equivalent, one has: 



$o (p) = <5o GVl+P 2 /(4m 2 )) • (4.19) 

Here we made use of the relation (3.6) between the relativistic and nonrelativistic momenta. We have 
therefore: 



N$ = lim 



2i(p - fc d ) 



p-,~k A 1 - itan [5 (p^l+p 2 /(^m 2 ))] 



2i I ^2m 2 (^/T+p 2 /m 2 - 1) - W2m 2 (^/l + k 2 /m 2 - 1) 



= lim 

p-^d 1 — itan [5o(p)\ 

+ (4.20) 

Here we have again used the relation (3.6) between the relativistic and nonrelativistic momenta. The 
result in the third line of eq. (4.20) is valid up to corrections of the order E^/m 2 . It can be rewritten 
in terms of the asymptotic S-wave normalization as follows: 



A S = A S 



1 + (1 + vi) 



2\ 



16m 



(4.21) 



To end this section, we note that the other deuteron properties such as in eq. (4.6), Qd in eq. (4.7)# 13 
and \J (r 2 )^ in eq. (4.8), which are not observable, are expected to change when calculated using the 
nonrelativistic wave function from eq. (4.3). For a recent reviews on the deuteron the reader is referred 
to refs. [70, 71]. 

* 13 It would be more appropriate to introduce a special notation for the quadrupole moment defined in eq. (4.7) using 
nonrelativistic impulse approximation in a way similar to the deuteron rms-radius. Unfortunately, no such notation 
appears in the literature, which might lead to a confusion. It should be understood that while the deuteron quadrupole 
moment represents the response of the deuteron to an external electromagnetic field and is certainly measurable, Eq. (4.7) 
gives only an approximation, which is model-dependent and not observable. 
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5 The fits 



In this section we discuss the determination and specify the values of the various LECs adopted in 
the present analysis. Throughout this work, we use the following values for the pion decay constant 
F w , the pion masses M n ±, M n o and the proton and neutron masses m p and m n : F w = 92.4 MeV, 
M n ± = 139.5702 MeV, M w o = 134.9766 MeV, m p = 938.2720 MeV, m n = 939.5653 MeV. 

For the 1PE potential, we use the expression (2.13) with gA = 1.29. This larger value of the LEC 
qa as compared to the standard one gA = 1-26 is in order to account for the Goldberger-Treiman 
discrepancy as discussed in section 2.3. Notice further that we take into account the pion mass 
difference as given in eq. (2.46). The leading 2PE potential given in eq. (2.14) is parameter-free. The 
NNLO and N 3 LO 2PE contributions in eqs. (2.15), (2.18)-(2.20) depend on the LECs c ± , c 2 , c 3 and c 4 
from the second-order ttN Lagrangian as well as on d± + d 2 , d 3 , and d\± — d\^ from the third-order 
ttN Lagrangian. For the LECs we adopt the central values from the Q 3 -analysis of the ttN system 
[73]: ci = -0.81 GeV -1 , c 4 = 3.40 GeV" 1 . For the constant c 3 the value c 3 = -3.40 GeV -1 is used, 
which is on the lower side but still consistent with the results from reference [73]: c 3 = —4.69 ± 1.34 
GeV" 1 . The same value for c 3 has been adopted in our NNLO analysis [21]. Further, this value was 
found in Ref. [26] to be consistent with empirical NN phase shifts as well as the results from dispersion 
and conventional meson theories. Notice, however, that it is about 25% smaller in magnitude than 
the value extracted from the partial wave analysis of the pp and np data [72]. The LEC c 2 could 
not be fixed accurately analyzing pion-nucleon scattering inside the Mandelstam triangle in [73]. We 
therefore adopt the central value found in the third-order analysis [23]: c 2 = 3.28 GeV -1 . For the 
combinations of d^s, we again use the values found in [23]: d-y + d 2 = 3.06 GeV" 2 , d 3 = -3.27 GeV" 2 , 
d 5 = 0.45 GeV" 2 and d 14 - d 15 = -5.65 GeV~ 2 . 

We now turn to short-range contact interactions. The two v = -LECs Cs,t, seven v = 2 -LECs 
C1...7 as well as fifteen v = 4 -LECs D1...15 in eq. (2.4) are unknown and have to be fixed from a fit 
to data (i.e. to Nijmegen phase shifts). Contributions of the contact interactions to various partial 
waves are given in eq. (2.5). Thus, we have to determine 8 LECs in the 3 S± — 3 D\ channel, 4 LECs 
in the 1 S channel, 3 LECs in the 3 P 2 - 3 F 2 channel, 2 LECs in each of the 1 P±, 3 Pi, 3 P partial 
waves and 1 LEC in each of the 1 D 2 , 3 D 2 , 3 D 3 — 3 G 3 channels. In addition to the above mentioned 
isospin-conserving contact interactions, we have two isospin-violating contact terms with unknown 
coefficients, see eqs. (2.35) and (2.37). Both terms contribute to the 1 So partial wave and provide 
charge-dependent contributions to the LEC Ciso- In the following, we will therefore distinguish 
between Cfg , C™^ and C™^ . We also note that we always use the proper kinematics as given in 
appendix A. Let us now specify precisely our way of fixing the LECs. The LECs contributing to 
isovector channels 3 P 2 — 3 F 2 , 3 Pi, 3 Po an d l D 2 have been fixed from a fit to Nijmegen pp phase shifts 
[38], which are much more precise than the corresponding np phase shifts. The isovector np phase 
shifts are then extracted from the pp ones in a parameter-free way by taking into account the proper 
1PE potential and switching off the electromagnetic interaction. This is precisely the same procedure 
as used in the Nijmegen PWA [38]. In the 1 Sq partial wave we have to take into account isospin- 
violating contact interactions as discussed above. We determine the LECs Cfg Q , C\so, D\so 
and D^ so from a combined fit in the 1 So pp and np channels. The LEC C 1 ™^ is then obtained from 
the requirement to reproduce the experimental value [74, 75] a nn = —18.9 fm for the nn scattering 
length. All remaining LECs are fixed from a fit to np phases from Nijmegen PWA [38]. We notice 
that contrary to our NLO and NNLO analysis [16, 21], we had to use here a large energy interval, i.e. 
up to Piab = 200 MeV, in order to fix the LECs. This is because of two reasons: first, the phase shifts 
in the 1 So and in the 3 5i- 3 I?i channels simply do not show enough structure beyond Ej a b = 100 MeV 
in order to fix reliably 4 and 8 parameters, respectively. Secondly, phase shifts at low energy are not 
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very sensitive to higher-order contact interactions except maybe in the two S-waves. 

It remains to specify the values for the cut-offs A and A which enter the Lippmann-Schwinger equation 
and the spectral-function representation of the two-pion exchange potential, respectively. Certainly, 
both cut-offs are introduced in order to remove high-momentum components of the interacting nucleon 
and pion fields, which are beyond the range of applicability of the chiral EFT. We remind the reader 
that from the formal point of view, one can choose any value for the SFR cut-off which is large 
enough so that the relevant physics is still present. Even the choice A = oo, which is equivalent to 
dimensional regularization, is formally possible since all terms with positive powers of A (and oc In A) 
can be absorbed by redefinition of the corresponding LECs. It has been argued in [21], however, that 
the choice A = 500 . . . 700 MeV leads to a natural separation of the long- and short-range parts of the 
nuclear force and allows to improve the convergence of the low-momentum expansion. In the present 
analysis we use this range for A. 

While A is related to perturbative renormalization of the pion loop integrals, the cut-off A specifies the 
way of nonperturbative renormalization of the Lippmann-Schwinger equation. Contrary to the SFR 
cut-off A, one, in general, cannot arbitrarily increase the value of A [76, 77, 78, 79]. This is because one 
needs an infinite number of counter terms in order to absorb all divergences arising through iteration 
of the potential in the Lippmann-Schwinger equation.* 14 Keeping A finite and of the order of the 
separation scale in the NN problem, one expects the contribution of the higher-order counter terms to 
be suppressed by powers of the generic low-momentum scale provided that the corresponding LECs 
are of natural size [76, 77]. In our previous NLO and NNLO analyses [81] based on the dimensionally 
regularized expressions for the potential we have used A = 500 . . . 600 MeV with the regulator function 
being defined as f A (p) = exp[— p 4 /A 4 ]. In the more recent study [21] based on the SFR approach we 
have increased this range to A = 450 . . . 650 MeV using f A (p) = exp[— p 6 /A 6 ]. We have, however, 
found in [21], that the upper values of A are already rather close to its critical value A c , above which 
one encounters spurious deeply-bound states. Notice that the values of various LECs start to strongly 
vary for A ~ A c leaving the natural range. In order to avoid such a situation we slightly reduce the 
range of variation of A to 450 . . . 600 MeV in the present analysis. To be more specific, we will use the 
following cut-off combinations (all values in MeV) : 

{A, A} = {450, 500}, {600, 600}, {450, 700}, {600, 700} . (5.1) 

For A = 500 MeV the value A = 600 MeV is already found to be close to A c . We therefore replace the 
cut-off combination {600, 500} by {600, 600}. Notice that further reducing of the A-values beyond 
A = 450 MeV is, in principle, possible but leads to a strong increase of the theoretical uncertainty. We 
therefore refrain from doing that. Finally, we notice that a more elegant regularization prescriptions, 
like e.g. lattice regularization, would allow to regularize pion loop integrals and the Lippmann- 
Schwinger equation in the same way without introducing two independent scales A and A. For a 
related recent discussion of the role and optimal choice of the cut-off A in the LS equation the reader 
is referred to Refs. [77, 79]. 

Let us now give the precise definition of the phase shifts considered in the present work and remind 
the reader on the type of phase shifts used in the Nijmegen PWA [38]. We will adopt here the notation 
of ref. [50] and denote by 5^ the phase shift generated by the potential W with respect to the solution 
with V as the interaction. 

* 14 It has been shown in [80] that l/r n singular potentials, which arise e.g. from pion exchange contributions, can be 
renormalized by a one-parameter square-well counterterm, see [8] for a related work. Although the authors of [80] have 
demonstrated that the low-energy NN observables can be made independent of the square-well width by adjusting the 
square- well strength, the power counting scheme adopted in the present work is not consistent with such an approach. 
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• pp phases. 

The full phase shifts <5em+n °f electromagnetic plus strong interaction can be expressed as: 

<5em+n = <5em + fi+N • (5-2) 

The above expression applies to uncoupled channels. For coupled channels one has to translate 
the addition law for the phase shifts into a multiplication law for the corresponding S-matrices 

•S'em+n = (<Sem) 1 ^ 2 <5em+n('S'em) 1//2 , (5.3) 

see [51, 82] for further details. Such a modification for coupled channels will, however, not change 
the conclusions of this section. We will therefore not consider the coupled case in detail in what 
follows. The last term in eq. (5.2) corresponds to the phase shifts of the electromagnetic plus 
nuclear interaction with respect to electromagnetic wave functions. These phase shifts ^f^+N 
are the ones which are given in the Nijmegen PWA [38]. Notice that the electromagnetic phase 
shift Sem can be represented as 

X X l XC1 | cCl+C2 . rCl+C2+MM c , / /r A\ 

OEM = 0C1 + C1+C2 + <3 C 1+C2+MM + ^C1+C2+MM+VP = d Cl + P + <P + T , (5.4) 

Cl I C2 Cl I C2 I MM 

where we have introduced the abbreviations p = c>ci+C2> <j> = ^ci+C2+mm ano - r = ^ci+C2+mm+vp - 
In practice, the quantities p, <\> and r are usually calculated using the CDWBA. This is justified 
due to the smallness of the corresponding interactions Vc2, Vmm and Vyp. In this case one has 
approximately: 

4> ~ ^Cl+MM , T ~ ^Cl+VP • (5-5) 

For more details the reader is referred to refs. [50, 51]. 

The phase shifts <5em+n obtained in the Nijmegen PWA do, however, not correspond to the type 
of phase shifts, which is usually considered in practical calculations, namely the phase shifts 
SqI , N of the modified Coulomb plus strong interactions with respect to the phase shifts of the 
modified Coulomb potential. These phase shifts can easily be calculated for any given nuclear 
potential using e.g. the method described in appendix C. The pp phases considered in the present 
work are of that type. We therefore need to relate ^ci+n ^° ^ ne previously discussed phase shifts 
^em+n- This can be done by noting that the total phase shift 5em+n can be expressed in the 
form: 

<>EM+N = <$C1 + <^C1+C2+MM+VP+N ■ ( 5 - 6 ) 

In the coupled case one has to modify this relation in a way analogous to eq. (5.3). Again, due 
to the smallness of the potentials Vq2, Vmm and Vyp one can make use of the DWBA to relate 

<5 C1+C2+MM+VP+N and 5 C1+N> which leads to [ 50 ] : 

(^ci+C2+mm+vp+n)« ~ (<5ci+n)/ = (5.7) 

= (^C1+N+C2+MM+Vp)' 

drxi{r) [V C 2{r) + Vuu{r) + Vyp(r)] Xl{r) , 



— m 



v 



where I is the angular momentum and Xi{ r ) 1S the wave function for the potential Vci + V^. 
Combining now eqs. (5.2)— (5.7) we end up with the following formula which relates the phase 
shifts <5q{ +n we are calculating to the ones 5™+N of the Nijmegen PWA: 



(SgJ+N)* ~ -h + n + h + Ti. (5.8) 
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This formula can be further simplified if one notes that for all Z > 1 the wave functions Xi( r ) near 
the threshold are almost not affected by the nuclear interaction [50]. Thus, one can approximately 
replace in eq. (5.7) the wave functions Xi{ r ) by the regular Coulomb wave functions, which leads 
to A; ~ pi + 4>i + 77. Therefore, for all partial waves except 1 S'o one has: (5q} +n )z ~ (^em+n)^- 
In the case of the 1 Sq partial wave, it has been argued based on explicit calculations that the 
quantity Ao is "sufficiently model-independent" for a wide range of nuclear forces at least in the 
case Vc2 = 0. Notice that the magnetic moment interaction in eq. (2.55) does not contribute in 
that and all other spin-singlet channels. The values of Ao based on the Nijmegen N78 potential 
[83] as well as for the tq and po are given in [50] at various energies. In the present analysis, we 
will use these values for the above mentioned quantities in order to relate our 1 Sq phase shift to 
the one of the Nijmegen PWA via eq. (5.8).* 15 Notice that a more accurate way to determine 
the phase shifts (<5§{ +N )i would be to define the regularized expressions for electromagnetic 
interactions and to calculate the quantities p, </>, r and A; explicitly (one can still use DWBA). 
We, however, believe that there is no need for such a refinement at the level of accuracy of N 3 LO. 

Let us now summarize our way of calculating phase shifts in the pp system. We compute the 
phase shifts <5ci+N of the modified Coulomb plus nuclear potential Vci + Vn with respect to wave 
functions of the Vci-potential. The strong interaction part of the chiral nuclear force at order 
N 3 LO is discussed in section 2. It has to be used in the relativistic Schrodinger equation (A. 4), 
while the electromagnetic interactions in eqs.(2.52), including the modified Coulomb force Vci, 
are to be used in the nonrelativistic-like equation (3.12). We therefore first apply eq. (3.15) 
to derive the modified strong potential for use in eq. (3.12).# 16 We then calculate phase shifts 
<^Ci+N m momentum space as described in appendix C. In the fitting procedure the calculated 
phase shifts ^ci+n are compared with the phases <5f{vi+N of the Nijmegen PWA [38]. For all 
partial waves except we use the approximation (5q} +n ); ~ (<^em+n)'- For * ne phase, we 
make use of equation (5.8), where the quantities A;, pi and t\ are taken from reference [50]. 

• np phases. 

As already pointed out in section 2.4.2, the electromagnetic interaction in the np case is given 
entirely in terms of magnetic moment interaction Vmm(^p)- Consequently, eq. (5.2) takes the 
form 

(5.9) 

The np phase shifts of the Nijmegen PWA as well as in our analysis correspond to the phase 
shifts <5mm+n of nuclear plus magnetic moment interactions with respect to magnetic moment 
interaction wave functions. Notice that the term in Vmm(^p) in eq. (2.55) proportional to A gives 
rise to the so-called "class IV" isospin-breaking force [84], which mixes spin-singlet and spin- 
triplet states. The contribution of this term is very small and usually only taken into account 
when constructing the magnetic moment scattering amplitude, see e.g. [85]. In our analysis we 
make use of the standard approximation [85] 

«+n ~ *n, (5.10) 



# 15 The quantity A; in [50] does not contain the contribution due to magnetic moment interaction, which has been 
neglected in that work. In our analysis we only need to know A; explicitly for I — and therefore can still use the result 
of [50]. 

* 16 Notice that in principle, the modified potential should contain contributions due to electromagnetic interactions. 
Equation (3.15) should actually be applied to the sum of the strong interaction and electromagnetic potentials Vn 
and Vem- The modified potential would then contain pieces oc Vn, Km, (1/m) VnVn, (1/m) VnKiii and (1/m) 14m Km- 
Applying eq. (3.15) only to the strong potential Vn and adding the appropriate (modified) electromagnetic contributions, 
we thus miss terms oc (1/m) VixrVim. Such contributions are suppressed by a factor Q/m compared to the ^-exchange 
in eq. (2.51) and beyond the accuracy of the present calculation. 
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for all / / 0-states. Therefore and because of the fact that the magnetic moment interaction 
does not contribute to the 1 S'o channel, we have to take into account Vmm(^p) explicitly only in 
the 3 Si- 3 -Di partial wave. In that case the phase shift <^mm+n * s calculated by subtracting 5mm 
from (5mm+n, where the phase shifts 5mm are obtained using the Born approximation. 

• nn phases. 

As in the previously considered case of neutron-proton scattering, in the nn system one has 
to take into account only magnetic moment interaction. Decomposing the phase shifts as in 
eq. (5.9) one can make use of the approximation <J^m+n ~ in ah nn partial waves. As already 
stated before, this approximation is accurate for partial waves with / / 0. In the case of the 
1 Sq partial wave, the phase shift is still given by 5n since the long-range magnetic moment 
interaction does not contribute to this channel. Therefore, nn phase shifts in all partial waves 
correspond to 5n- 



6 Results and discussion 
6.1 Phase shifts 

In the following sections, we will show the results for the np phase shifts. Before showing the 
results of our analysis, let us make a simple estimate for the expected theoretical uncertainty at 
N 3 LO. Following the reasoning of ref. [21], we expect for the uncertainty of a scattering observable 
at c.m.s. momentum k at N 3 LO to be of the order ~ (max[fe, M T ]/X) 5 . To provide a fair esti- 
mate, we identify the hard scale A with the smallest value of the ultraviolet cut-off, i.e. we adopt 
A ~ 450 MeV. This results in the following estimations for the theoretical (maximal) uncertainty: 

~ 0.5% at E iah ~ 50 MeV and below, 

~ 7% at E lah ~ 150 MeV, 

~ 25% at E lab ~ 250 MeV. 

One should keep in mind that the above estimations are fairly rough. For a detailed discussion on the 
theoretical uncertainty, especially at NLO and NNLO, the reader is referred to [21]. 

6.1.1 S waves 

The phase shifts in the 1 Sq and 3 Si partial waves are shown in Fig. 3. Both are visibly improved 
compared to the NNLO result. For E\ a i, = 50 MeV, 150 MeV and 250 MeV we find the phase shift in 
the partial wave in the ranges 40.42° . . . 40.72°, 16.04° . . . 17.03° and 2.22° . . . 4.76°, respectively. 
These values agree well with the ones from the Nijmegen PSA: 5 = 40.54°, 5 = 16.94° and 5 = 1.96°. 
The relative uncertainty of our results is in agreement with the above estimations except for E\ a ^ = 250 
MeV, where the phase shift is close to 0. The results for the 3 5i partial wave are similar to the ones 
in the 1 So channel. The uncertainty due to the cut-off variation is found to be smaller in this case. 

The pertinent S-wave LECs are tabulated in table 2 for the four pairs of cut-offs (5.1). Here, several 
remarks are in order. First, we note that, in general, one has to expect multiple solutions for the 
LECs. This problem has already been discussed in [16] at NLO and NNLO. For the 1 So channel, we 
have to fix five LECs Cfgn, C™E , CiS0j E>\ so , Df SQ , from a fit to Nijmegen pp and np phase shifts. We 
did find multiple solutions for LECs which describe the data equally well if we neglect isospin breaking 
and fix C™g Q , C\so, D\ so , Df so from a fit to the Nijmegen np phase shift. Taking into account isospin 
breaking effects and performing a combined fit to both pp and np phase shifts turns out to improve 
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the situation and help to sort out the true solution. We found a single solution for the LECs for the 
cut-off combinations {A, A} = {450, 700} and {600, 700}. For the two other cut-off combinations 
several local minima in the x 2_ pl°t have been observed. We then adopted the values for the LECs 
corresponding to the global minimum. We have checked that these values result from the ones for 
different {A, A} by a continuous change of the cut-offs. 

The situation in the 3 Si— 3 -Di— channel is even more complex since one has to determine eight LECs. 
Our results for the LECs are shown in tables 2, 3 and 4. Due to the large dimension of the parameter 
space, we cannot definitely claim that the found values for the LECs correspond to a true global 
minimum of the x 2 - 

Let us now comment on the naturalness of the determined LECs. In general, the natural size for the 
LECs can be (roughly) estimated as follows: 
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(6.1) 



where Alec is the scale entering the values of the LECs and the factor Air results from the angular 
integration in the partial wave decomposition, see appendix B for the details. The S-wave LECs 
shown in table 2 are of the natural size except the LECs D\ SQ and D^ S1 , which are somewhat large 
in magnitude. Indeed, estimating the scale Alec as Alec ~ 500 MeV leads to \D{\ ~ 2.4 in the 
same units as used in table 2. Still, the higher-order contact interactions are suppressed compared to 
the lower-order operators at low momenta. For example, for the cut-off combination {450, 500} and 
p = p' = M-k the contributions of the contact operators at various orders are given by: 



p=p'=M 7r 



CTso + C\so(p 2 +P' 2 ) + (D\soP 2 p' 2 + D 2 1S0 ( P 4 +p> 4 : 



- 0.091 + 0.057 + (-0.010 + 0.003) 



x 10 4 GeV~ 2 . 



p=p'=M 7r 

(6.2) 



6.1.2 P waves 

Our results for the np P-waves and the mixing angle t\ are shown in Fig. 4. All phase shifts are 
visibly improved compared to the NLO and NNLO results. One has, however, to keep in mind that 
two independent parameters appear now in each of these channels instead of one parameter at both 
NLO and NNLO. The results of the Nijmegen PWA are reproduced in our N 3 LO analysis within the 
theoretical uncertainty in all phase shifts with exception of 1 P\ at larger energies. In the case of the 
3 P2 partial wave, the band is dominated by variation of the SFR cut-off. In particular, lower values 
for this cut-off lead to a better agreement with the data at higher energies. At £i a b = 250 MeV, the 
N 3 LO phase shifts deviate from the data by an amount of up to ~ 8°. The typical size of the P-wave 
phase shifts at this energy is of the order ~ 25° . The uncertainty in the calculated phase shifts due 
to the cut-off variation agrees therefore fairly well with the estimation in section 6.1. We remind the 
reader that the theoretical bands at NLO and NNLO are expected to have a similar width, since the 
effective potential at these orders contains the same set of contact interactions (counter terms). As 
explained in detail in [21], the uncertainty resulting from the cut-off variation at NLO is smaller than 
the actual theoretical uncertainty at this order. 

The pertinent P-wave LECs are tabulated in table 3 four these four pairs of cut-offs. As in the 
previously discussed case of the S-waves, we found multiple solutions for the LECs. The physical 
solution can easily be determined due to the smaller number of parameters (two unknown LECs in each 
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P-wave). To illustrate this point consider the 3 Pi partial wave with the cut-offs {A, A} = {450, 500}. 
We find two solutions for the LECs C 3P \ and D 3P i fitting to the Nijmegen pp phase shifts: 



where we used the same units as in Table 3. Both sets of parameters lead to an accurate description 
of the data, which is shown in Fig. 5. The solution in the first line of eq. (6.3) satisfies the naturalness 
assumption for the LECs and has been adopted in the present analysis. Notice further that the value 
C 3P1 = -0.6334 x 10 4 GeV -2 is close to the NLO and NNLO values for these LECs (for the same 
cut-off combination) C 3P1 = -0.4932 x 10 4 GeV~ 2 and C 3P1 = -0.7234 x 10 4 GeV^ 2 , respectively. 
The results for other partial waves are similar. All LECs Cj in the P-waves are found to be of natural 
size and take the values which are close to the ones at NLO and NNLO. The P-wave LECs Di are 
natural as well. 

Let us now comment on isospin-breaking. As already explained before, the LECs in the isovector 
partial waves are fitted to the pp phase shifts. To calculate the corresponding np phase shifts, we 
switch off the electromagnetic interaction and adjust for the proper pion mass in the 1PE potential. 
The differences between the corresponding pp and np phase shifts at three different energies E^ = 10, 
25 and 50 MeV are shown in Tables 5, 6 and 7, respectively. In general, we see that the effects 
due to removal of the Coulomb interaction, Aj, agree very well with the ones of Nijmegen PWA. 
The uncertainty due to the cut-offs variation becomes larger at higher energies. The effects due to 
including the pion mass difference, AJ , show typically somewhat larger deviations from the Nijmegen 
PWA. This is presumably to a large extent due to a different treatment of the 1PE force: while 
we use the potential in momentum space with high momenta being cut off, the Nijmegen group 
performs calculations in coordinate space, and chooses to cut-off the long-range potential at 1.4 fm. 
Although both methods certainly lead to the same long-distance asymptotics of the 1PE potential, 
they differ significantly in the treatment of its shorter-range part. The largest deviations from ~ 15% 
at £ lab = 10 MeV to ~ 20%-30% at E Xab = 50 MeV for Af from the Nijmegen PWA are observed 
in the 3 P2 partial wave. It is comforting to see that both isospin-violating effects (i.e. due to the 
Coulomb force and the pion mass difference in 1PE) are in most cases of the same size, as it is also 
expected from power counting arguments, see section 2.4 for more details. We also note that these 
two effects have often opposite sign and tend to cancel. For example, one observes at E^b = 10 MeV 
[38]: A 3P0 = A^ po + AJ P0 = 0.371° - 0.447° = -0.076°. 

6.1.3 D— and higher partial waves 

The results for D-, F- and G-waves are shown in Figs. 6, 7 and 8, respectively. We remind the reader 
that at N 3 LO there is one adjustable constant Di in each of the D-waves, while F- and higher partial 
waves are parameter-free. As depicted in Fig. 6, the shape of the 3 D 3 partial wave is still not properly 
reproduced at N 3 LO, although it is greatly improved compared to NLO and NNLO predictions. This 
phase shift is, however, rather small as compared to other D-wave phase shifts, and thus one expects 
relatively small effect of this phase shift on the NN scattering observables. We also note that the 
absolute deviation from the data in this channel is not larger than in the other D-waves. Most of the 
F- and G-waves are at N 3 LO in agreement with the data. One observes that the theoretical bands 
do not get thinner at N 3 LO, which might at first sight appear strange. This, however, is naturally 
explained by the fact that there are no short-range contact terms in these channels. Such terms 
start to contribute to F-waves at N 5 LO (Q G ) and to G-waves at N 7 LO (Q 8 ). Consequently, one 
should expect the uncertainty due to the cut-off variation to be of the same size for calculations up to 



and 



C 3P1 = -0.6334 
C 3P1 = 5.9620 , 



D 3Pl = 4.2359 , 
D 3P1 = -20.6154, 



(6.3) 
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these high orders in the chiral expansion. Clearly, peripheral partial waves are strongly dominated by 
the 1PE potential, which represents the longest-range part of the strong nuclear force. Indeed, one 
observes that the phase shifts are mostly well reproduced already at NLO, while NNLO and N 3 LO 
corrections only produce minor changes. Notice further that due to the smallness of the phase shifts, 
the Born approximation works very well in high partial waves and the phase shifts are essentially given 
by the diagonal (in momentum space) matrix elements of the two-nucleon potential. It is then clear 
that the bands arise almost completely due to multiplying the potential by the regulator function. 
The only exception from this rule is given by the 1 Z?2 partial wave, where the SFR cut-off has a larger 
impact on the phase shift at higher energy than the cut-off in the Lippmann-Schwinger equation. In 
particular, the lower values of this cut-off lead to larger values of the phase shift. 

The determined LECs Dj are tabulated in Table 4. All of them are of natural size. 
6.2 S-wave effective range expansion 

We now regard the S-wave effective range parameters and begin with the np system. In that case 
one can make use of the usual effective range expansion for finite-range potentials, eq. (D.l). The 
reason is that the long-range magnetic moment interaction does not contribute to states with 1 = 0. 
Notice, however, that one should not use the standard effective range expansion for the 3 D\ partial 
wave and mixing angle e±, which are modified in the presence of the long-range (~ 1/r 3 ) magnetic 
moment interaction. Our results for the 3 <Si and 1 Sq scattering length, effective range and shape 
coefficients ^2,3,4 are summarized in Tables 8 and 9. The results for the 1 Sq scattering length and 
effective range are improved compared to the NLO and NNLO predictions of [21]. The N 3 LO result 
for the scattering length fills a small gap between the NNLO prediction and the value of the Nijmegen 
PWA. The uncertainty due to the cut-off variation for all effective range parameters turns out to be 
smaller at N 3 LO compared to NNLO, as it should. We observe a minor discrepancy for the shape 
coefficient i>2, which might however simply reflect the lack of numerical accuracy, with which this 
quantity is calculated. The description of the 3 Si effective range parameters is similar to the one in 
the 1 Sq channel. 

Next, we consider the pp system. This case is much more complex since one has to account for 
electromagnetic interaction. Ideally, one should use the phase shifts ^f^J+N anc ^ the expression for 
the effective range function given in [50] to obtain the effective range expansion for the nuclear force 
in presence of the long-range electromagnetic interactions, which in the 1 Sq channel are given by 
the improved Coulomb and vacuum polarization potentials in eqs.(2.53) and (2.56). In the present 
analysis we have used a simplified description for the pp phase shift as explained in section 5. We 
do not calculate explicitly the phase shifts <5™ +N but rather the ones c>c 1+ n °f nuclear plus modified 
Coulomb potential with respect to Coulomb wave functions, adjusting for the difference as explained 
in section 5. We have therefore made use of eq. (D.5) to calculate the pp scattering length and effective 
range. We obtain the following values: 

a vv = -7.795...- 7.812 fm, r pp = 2.73 ... 2.76 fm , (6.4) 

where the uncertainty is due to the cut-off variation. These values agree nicely with the experimental 
ones [82]: 

a c p f = -7.8149 ± 0.0029 fm , r^ p = 2.769 ± 0.014 fm . (6.5) 

One should, however, keep in mind that we made an approximation and neglected the effects due 
to the long-range part of the vacuum polarization potential and the interaction in the second line 
of eq. (2.53). As found in [50] neglecting these electromagnetic interactions affects the values of a pp 



32 



and r pp by an amount smaller than 0.01 fm, which is within the theoretical uncertainty of the present 
analysis. 

Finally, we consider the nn system. Since no long-range electromagnetic interactions contribute to the 
1 Sq partial wave, one can use the effective range expansion (D.l). Since we have used the "standard 
value" for the nn scattering length ci nn 

o^ = -18.9 ±0.4 fin (6.6) 

as an input to fix the LEC of the leading isospin-violating short-range interaction, we can only make 
predictions for the effective range r nn : 

r nn = 2.76... 2.80 fm. (6.7) 

This agrees with the experimental number [86] 

= 2.75 ±0.11 fm. (6.8) 

Notice that there is still some controversy about the experimental value of the nn scattering length 
extracted using different reactions. For example, the value a^r? = —18.50 ±0.53 fm has been reported 
from studying the 2 H(-7r~, ny) n process [74], while measurements of the neutron deuteron breakup 
reaction lead on one hand to ann = —18.7 ± 0.6 fm [87] and on the other hand to different values, 
a^v = -16.1 ± 0.4 fm and cwf = -16.3 ± 0.4 fm [88]. 

Last but not least, we would like to point out that the relation between the values of the effective range 
parameters, which would result if there would be no electromagnetic interaction, and the observed 
ones is highly nontrivial. Neglecting electromagnetic nucleon mass shifts, we can, for example, switch 
off the Coulomb interaction in the pp system and recalculate the effective range coefficients using 
eq. (D.l). This leads to: 

a pp = -16.00 ... - 16.63 fm , r pp = 2.81 . . . 2.86 fm . (6.9) 

Although the value for a pp is fairly close to the one for the nn and np scattering lengths, as one would 
expect from the approximate isospin invariance of the strong interaction, it should be understood that 
effects due to electromagnetic interaction are not completely removed from these quantities. To clarify 
this point let us take a look at the leading nn, np and pp short-range interactions: 



150 — V^lSOMr + Pn 



e 2 



(4tt) 2 ' 

~np ~n e 2 

^lSO = (ClfoW + Pnp ' 



e 2 



Cf P so = (Cr S0 Ur + P PP j^. (6.10) 

Here the LECs (Cl so ) s tr are entirely due to the strong interaction. If only linear terms in the quark 
mass difference are included, see eq. (2.35), these LECs are related with each other as (Cfg ) s tr + 
(Ciso)str = 2(C'™^ ) s tr and the difference (C , ^ ) str — (C'ljoW is proportional to eM 2 . The terms oc (3i 
in eq. (6.10) are due to the short-range electromagnetic interactions, see eq. (2.37).* 17 Since we do not 

# 17 The fact that we have three independent LECs f3 nn , [3 np and [3 PP and only two terms in eq. (2.37) might appear 
confusing. In fact, we have only shown explicitly electromagnetic isospin-breaking and omitted isospin conserving terms 
in eq. (2.37). One of the two electromagnetic isospin conserving contact interactions contribute to 1 So and one to 3 Si 
NN scattering. Therefore, three and not two independent electromagnetic terms contribute to 1 So NN scattering. 
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know the values of the LECs fa in eq. (6.10) and it is not possible to disentangle them from C{ so in the 
two-nucleon system, we cannot extract the values for NN observables due to the strong interaction out 
of the experimentally measured quantities. Notice that the LECs /3j might (at least in principle) be 
determined from processes with external pions. Notice further that the pp scattering length with the 
long-range Coulomb interaction being switched off is even not a well-defined quantity in an effective 
field theory approach since it is sensitive to details of the strong interaction at short distances. Indeed, 
the extracted scattering length a pp in eq. (6.9) shows a significant cut-off dependence. Clearly, the 
scattering length due to pure strong interaction is perfectly well defined and the cut-off dependence is 
(largely) absorbed by the appropriate "running" of [3 PP . For related discussion on the proton-proton 
scattering length in context of effective field theory see [89, 90]. Furthermore, a useful approximation 
for the quantity a pp based on the short-range nature of the strong interaction and the fact that the 
scattering length is large can be found in [91]. 

6.3 Two-nucleon scattering observables 

Once the NN phase shifts are calculated, all two-nucleon scattering observables can be obtained in 
a straightforward way using e.g. the formulae collected in [52]. In Figs. 9, 10, 11 and 12 we show 
the np differential cross section and vector analyzing power at E^ =25, 50, 96 and 143 MeV at 
NNLO and N 3 LO in comparison with the data and the Nijmegen PWA results. In this calculation, 
we have included all np partial waves up to j < 8 and did not take into account the magnetic moment 
interaction. At the lowest energy we have considered, E^ = 25 MeV, both NNLO and N 3 LO results 
are consistent with the ones of the Nijmegen PWA. The small disagreement with the Nijmegen PWA 
in the analyzing power at forward direction is due to the neglected magnetic moment interaction. At 
higher energies the NNLO predictions become less precise. At N 3 LO the uncertainty in the cross- 
section due to the cut-off variation at the largest energy we have calculated, E^b = 143 MeV, is 
less than 10%. It is comforting to see that NNLO and N 3 LO results overlap in most cases and are 
both in agreement with the Nijmegen PWA. We further notice that the small but visible deviations 
of our N 3 LO result for the differential cross section from the Nijmegen PWA curve at forward and 
backward angles and higher energies is most probably due to the lack of partial waves with j > 8 in 
our calculations. The convergence of the partial wave expansion is well known to be slow in these 
particular cases. For example, it has been found in [92] that a sum up to j = 16 is needed to obtain 
convergence for the cross section at E^ = 300 MeV within 1%. 

6.4 Deuteron properties 

We now turn to the bound state properties. We stress that we do not use the deuteron binding energy 
as a fit parameter as it is frequently done but rather adopt the same LECs as obtained in the fit to 
the low phases. In Table 10 we collect the resulting deuteron properties in comparison to the NLO 
and NNLO results from [21]. All results for the deuteron properties in this table have been calculated 
using the formulae given in section 4 based on the relativistic wave function ^ d (p). First, we note a 
clear improvement at N 3 LO in the chiral expansion. The predicted binding energy at N 3 LO is within 
0.4% of the experimental value. This has to be compared with 1%-1.5% (~2%-2.5%) deviation at 
NNLO (NLO). Also visibly improved is the asymptotic S-wave normalization strength As, which 
now deviates from the experimental (central) value by 0.3% as compared to ~1.1% (~1.9%) at NNLO 
(NLO). Our predictions for the asymptotic D/S-ratio have a tendency to slightly reduce its value when 
going from NLO to NNLO to N 3 LO. The results at all orders are in agreement with the data within the 
experimental uncertainty. Further, our N 3 LO result for rj^ agrees well with the one of the Nijmegen 
PWA [93], r/d = 0.0253(2). We do not observe any improvement for the quadrupole momentum at 
N 3 LO, which shows an even larger deviation from the data compared to NNLO (6% - 8% versus 4% - 
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5%). We, however, remind the reader that the present calculation of Qd is based on formulae of section 
4. It is incomplete and does, in particular, not take into account the contribution of the two-nucleon 
current. Notice that apart from the pion-exchange two-nucleon currents, there are contributions from 
two-nucleon contact current, where the corresponding LEC cannot be fixed from nucleon-nucleon 
scattering. Such current results from the operator in the effective Lagrangian with four nucleon fields, 
one photon field and two derivatives. It appears natural to fix the value of the accompanying LEC from 
the requirement to reproduce the value of the deuteron quadrupole moment. For the calculations of the 
various deuteron properties including the quadrupole moment as well as other two-nucleon observables 
in pionless EFT the reader might consult Refs. [94, 95]. Notice that the situation with the quadrupole 
moment is analogous to the one described in [96] for the deuteron magnetic moment. In that case 
the corresponding short-range two-nucleon current results from the operator with just one derivative 
and thus appears even at a lower order. The situation with the deuteron rms-radius is similar to the 
one with the quadrupole moment: we observe a larger deviation from the data at N 3 LO as compared 
to the NLO and NNLO results. Notice however that the deviations at N 3 LO from the experimental 
number are still of the order of 0.5% or less. The above comment on the missing contributions in 
the quadrupole moment calculation applies to the deuteron rms-radius as well. At N 3 LO one should 
account for the contribution due to the short-range two-nucleon current, which results from the contact 
operator with four nucleon fields, one photon field and two derivatives, see [94] for more details. It is 
remarkable that the cut-off dependence of \J (r 2 )^ at N 3 LO is significantly larger compared to NLO 
and NNLO. This implicitly confirms our previous statement about the necessity to incorporate the 
short-range current at this order. The cut-off dependence of the corresponding LEC will compensate 
the cut-off dependence of \J (r 2 )^ making the deuteron "point-nucleon" radius cut-off independent 
up to higher-order corrections. The complete N 3 LO calculation of the quadrupole moment and the 
"point-nucleon" electric charge radius of the deuteron will be presented in a separate publication. As 
a numerical check, we have recalculated all deuteron properties using the nonrelativistic wave function 
4> d (p) in eq. (4.3). As expected from the discussion in section 4, we reproduce the values for and 
r/d- The asymptotic S-wave normalization As changes by A$ — As = 0.00039 fm -1 / 2 . This has to be 
compared with the value As — As = 0.000392 fm -1 / 2 from eq. (4.21). The quadrupole moment and 
the rms-radius change by 0.4% and 0.1%, respectively. Finally, we show the deuteron wave function in 
coordinate space in fig. 13 for a particular cut-off choice, together with results obtained at NLO and 
NNLO. One observes a stronger suppression of the S-wave component at short distances compared to 
NLO and NNLO, as well as the lower probability for the deuteron to be in the D-state. The latter 
observation follows also from the smaller value of quoted in Table 10. The shape of the wave 
function changes for different cut-off choices. We remind the reader that the deuteron wave function 
is not observable (except at very large distances). 

7 Summary 

In this paper, we have considered the interactions between two nucleons at N 3 LO in chiral effective 
field theory. The pertinent results of this study can be summarized as follows: 

i) The two-nucleon potential at N 3 LO consists of one-, two- and three-pion exchanges and a set 
of contact interactions with zero, two and four derivatives, respectively, according to the chiral 
power counting, see also table 1. We have applied spectral function regularization to the multi- 
pion exchange contributions. This allows for a better separation of the low and high momentum 
components in the pion loop diagrams than dimensional regularization. Within this framework, 
we have shown that three-pion exchange can safely be neglected. The corresponding cut-off is 
varied from 500 to 700 MeV. The LECs related to the dimension two and three NNtttt vertices 
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are taken consistently from studies of pion-nucleon scattering in chiral perturbation theory, 
[23, 73]. In the isospin limit, there are 24 LECs related to four-nucleon interactions which feed 
into the S-, P- and D-waves and various mixing parameters, cf. Eq. (2.5). 

ii) We have reviewed the various isospin breaking mechanisms and proposed a novel ordering 
scheme, based on one small parameter that collects strong as well as electromagnetic isospin vio- 
lation, cf. Eq. (2.32) accompanied by a particular counting rule for photon loops, see Eq. (2.33). 
This differs from the scheme proposed and applied in Ref. [7]. In the actual calculations, we 
have included the leading charge-independence and charge-symmetry breaking four-nucleon op- 
erators, the pion mass difference in the 1PE, the kinematical effects due to the nucleon mass 
difference and the same electromagnetic corrections as done by the Nijmegen group (the static 
Coulomb potential and various corrections to it, magnetic moment interactions and vacuum 
polarization). This is done because we fit to the Nijmegen partial waves. In the future, it would 
be important to also include isospin violation in the 2PE, ^-exchange and the isospin breaking 
corrections to the pion-nucleon scattering amplitude (which have been consistently determined 
in [44]). 

iii) We have discussed in some detail the form of the scattering equation that is used to iterate 
the potential and similar for the bound state. We use the Lippmann-Schwinger equation with 
the relativistic form of the kinetic energy. Such an approach can easily be extended to external 
probes or few-nucleon systems. We have also discussed the reduction to a nonrelativistic form 
which be might of easier use in some applications. The LS equation is regulated in the standard 
way, cf. Eq. (2.28), with the cut-off varied from 450 to 600 MeV. 

iv) The total of 26 four-nucleon LECs has been determined by a combined fit to some np and pp 
phase shifts from the Nijmegen analysis together with the nn scattering length value a nn = 
— 18.9 fm, as detailed in section 5. The resulting LECs are of natural size except D\ SQ and D^ S1 . 
Comparing to the fits at NLO and NNLO, we had to extend the fit range to higher energies for 
the reasons discussed in section 5. 

v) The description of the low phase shifts (S, P, D) is excellent, see Figs. 3-6. In all cases, the 
N 3 LO result is better than the NNLO one with a sizeably reduced theoretical uncertainty. This 
holds in particular for the problematic 3 Po wave which was not well reproduced at NNLO. The 
peripheral waves (F, G, H, . . .), that are free of parameters, are also well described with the 
expected theoretical uncertainty related to the cut-off variations, see Figs. 7-8. We stress that 
the description of the phases in general improves when going from LO to NLO to NNLO to 
N 3 LO, as it is expected in a converging EFT. 

vi) The resulting S-wave scattering lengths and range parameters in the np (cf. tables 8 and 9) and 
pp systems (cf. Eq. (6.4)) are in good agreement with the ones obtained in the Nijmegen PWA. 
In addition, we can give theoretical uncertainties for all these quantities, which are mostly in 
the one percent range. 

vii) The scattering observables (differential cross sections, analyzing powers) for the np system dis- 
played in Figs. 9-12 are well described, with a small theoretical uncertainty at the order consid- 
ered here. 

viii) The deuteron properties are further predictions. In particular, we have not included the binding 
energy in the fits, the deviation from the experimental value is in the range from 0.4 to 0.07%. 
The asymptotic S-wave normalization and the asymptotic D/S are also well described. The 
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remaining discrepancies in the quadrupole moment and the rms matter radius are related to the 
short-ranged two-nucleon current not considered here. 

In the future, these studies should be extended in various directions. In particular, one should construct 
the electroweak current operators to the same accuracy and work out the corresponding three-nucleon 
force, which is of special interest since it does not contain any novel LECs. Furthermore, a more 
systematic study of isospin violation in the two- and three-nucleon systems based on the formalism 
developed here should be pursued. Work along these lines is under way. 
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A Kinematics 



Consider two nucleons moving with momenta p\ and p2 . We use the relativistic kinematics for relating 
the energy E^ of two nucleons in the laboratory system to the square of the nucleon momentum p 
in the center-of-mass system, which is defined by the condition p\ +P2 = 0.* 18 The relation between 
Slab and p 2 reads (here and in what follows: p = \p\): 



• Proton-proton case: 



Neutron-neutron case: 



Neutron-proton case: 



p 



P 2 = 7;m p E Uh . 



2 1 

p = -m n Ei ah . 



m 2 p E^ h {E^ h + 2m n ) 



(m n + m p ) 2 + 2E\ ab m p ' 

The relativistic Schrodinger equation for two protons or two neutrons in the cm. system reads: 

{2\Jp 2 + m - 2rnj + V * = E^ , 
where m is the proton or neutron mass. For the neutron-proton system it takes the form 



VT+m| + \Jp 2 +m 2 — m n — m p j + V 
The free Hamiltonian Hq can be expressed in terms of the mass m defined as 



m 



2m p m n 
m p + m n 



(A.l) 
(A.2) 

(A.3) 

(A.4) 
(A.5) 
(A.6) 



* 18 It would be more appropriate to call such a system center-of-momenta or rest-frame and not center-of-mass as 
usually done in the literature. 



37 



in the following way 



H = v / P 2 + m 2 + y /p 2 



+ ml-m n - m t 



2m , 



(A.7) 



modulo terms which are proportional to (m p — m n ) 2 . Taking into account such terms goes beyond 
the accuracy of the present analysis. We will therefore use the approximate expression (A.7) in this 
work, which leads to the Schrodinger equation of the type (A. 4). 



B Partial wave decomposition of the NN potential 

In this appendix we describe the partial wave decomposition of the two-nucleon potential. For that 
we first rewrite the potential V in the form 



V = V c + V a <?i • a 2 + V S l i 2(^1 + a 2 ) ■ (k x q) + V aL a x - (qx k) a 2 - (qx k) 

+ Vaq {#1 ■ 0) {#2 ■ 0) + V ak (<Tl ■ k) {B 2 ■ k) , 



(B.l) 



with six functions Vc(p,p', z), V a k{p,p' , z) depending on p = \p\, p' = \p'\ and the cosine of 

the angle between the two momenta is called z. These functions may depend on the isospin matrices 
r as well. To perform the partial wave decomposition of V, i. e. to express it in the standard Isj 
representation, we have followed the steps of ref. [97]. In particular, we start from the helicity state 
representation \p\i\ 2 ), where p = p/p and Ai and A2 are the helicity quantum numbers corresponding 
to nucleons 1 and 2, respectively. We then expressed the potential in the |jmAiA2) representation using 
the transformation matrix (p \i\ 2 \jm\i\ 2 ) , given in ref. [97]. The final step is to switch to the \lsj) 
representation. The corresponding transformation matrix (lsjm\jm\iX 2 ) is given in refs. [98], [97]. 

For j > 0, we obtain the following expressions for the non-vanishing matrix elements in the \lsj) 
representation: 



imv\m 
(mv\m 



2ir 



2ir 



I" dz [V C - W a +p' 2 p 2 (z 2 - l)V aL - q 2 V aq - k 2 V ak } P 3 (z) , 



V c + V a + 2p'pzV S L ~ P' 2 P 2 (1 + 3z 2 )V aL + 4k 2 V aq + -q 2 V ak 



x PAz) + 



-p'pVsl + 2p' 2 p 2 zV aL - 2p'p (V aq - ^V a; ) 



X {Pj-^Z) + P j+1 (z)) 



(j ± 1, lj\V\j ± 1, \j) = 2irJ^dz I p'p 



-V SL ± 



-p'pzV aL + V aq - -V c 



(j±i,ii|y|j T i,ij) = 



xPj(z) 



2 j + l V ' - aq i 

V C + V a +p'pzV SL +P ,2 P 2 (1 - z 2 )V aL 



ak 



Pj±i(z) 
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+ 



+ 



2j + 

±^ J ^V aL +p 2 (2V aq + l -V ak ) P j± i(*)|. (B.2) 




Here, Pj(z) are the conventional Legendre polynomials. For j = the two non-vanishing matrix 
elements are 

(000|V|000) = 2tt J' dz [v c - 3V a +p' 2 p 2 {z 2 - l)V aL - q 2 V aq - k 2 V ak } , 
(110|V|110> = 2vr J dz^zV c + zV a +p'p(z 2 -l)V SL +p' 2 p 2 z(l-z 2 )V aL 

- ((p' 2 + P 2 )z - 2p'p) V aq - (V 2 + p 2 )z + 2p'p) V ak } . (B.3) 

Note that sometimes another notation is used in which an additional overall minus sign enters the 
expressions for the off-diagonal matrix elements with I = j + 1, I' = j — 1 and I = j — 1, I' = j + 1. 

C Momentum space treatment of the Coulomb interaction 

In this appendix we would like to explain our way of treating the nucleon-nucleon scattering problem in 
the presence of the Coulomb interaction in momentum space (following closely Ref . [7] ) . The starting 
point is the nonrelativistic Schrddinger equation of the form (3.7) or (3.12), where the the potential 
consists of two pieces: the short-range one given by the strong interaction and the long-range one 
given by the Coulomb force. 

As the Coulomb potential is of infinite range, the S-matrix has to be formulated in terms of asymptotic 
Coulomb states. Therefore, the phase shifts for a given angular momentum I due to the strong potential 
in the presence of the long-range electromagnetic interactions, denoted by S\, are defined in terms of 
a linear combination of (ir)regular Coulomb-functions F(G) as 

xi(r)=F,(r)+tan(5i)G,(r) (C.l) 

analogously to the expression for an arbitrary potential of short range (i.e. in the absence of the 
Coulomb force) 

X?(r)=if(r)+tan(*f)G?(r) (C.2) 

with F°, G° denoting solutions of the Coulomb problem with zero charge (conventionally expressed 
in terms of Bessel and Neumann functions) and the corresponding phase shift is called 8? . So far, we 
have restricted ourselves to uncoupled channels. We will consider the coupled case later on. 

As eq.(C.l) exhibits asymptotical Coulomb-states, we have to re-express our Lippmann-Schwinger- 
equation in terms of them. A very convenient scheme for inclusion of the Coulomb force in momentum 
space was suggested long time ago by Vincent and Phatak [99] and is used in the present analysis. In 
what follows we will briefly describe this approach. 

The starting point of this technique is the observation, that for a potential of the form 

V = V c + V s (C.3) 
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with 

Vs^i = (r>R), (C.4) 

and ipi the two-nucleon wave-function for a given angular momentum, two exact solutions for the wave- 
function can be given for every point on a sphere with radius R + e. One is of the form as in eq.(C.l), 
and another one according to eq.(C.2) with the phase shifts calculated for the following potential as in 
eq.(C.3), with Vc, however, being the Fourier-transformed Coulomb-potential integrated to the radius 

R, 



V c (\q'-q\) = f (FreW-V*- 
Jo r 



t^ttoCI -cos(| q'-q\R)) . (C.5) 



Here, q, q' are the cms momenta and a is the fine-structure constant. On the above-defined sphere, 
both wave functions describe the same system. Now we know how to obtain an expression for the 
strong phase shift 8\ in the presence of the Coulomb interaction in terms of the short-range shift Sf in 
the absence of electromagnetism: We only have to match the two solutions. This is most conveniently 
done by requiring the logarithmic derivative of both solutions to be equal, what enables us to express 
the strong shift in the presence of the Coulomb force in a Wronskian form: 

, _ tan(^)[F,G ] + [F,F ] 

tan( ^ } " [F ,G]+tan(«5f)[G ,G] (C - 6) 

with 

™ = ( c £- p *L (c - 7) 

Let us now extend the previous consideration to the coupled case. For that we replace eq. (C.l) by 
the matrix equation: 

x l (r) = F(r) - mqK l G(r) , (C.8) 

where K l is the K-matrix for the strong potential in the presence of the Coulomb interaction, x l ( r ) 
is the 2x2 matrix which contains the wave functions 

y l (r) = ( xJ-^-it 7 ") x)-i,j+i( r ) \ , CQ) 
and F(r) and G(r) are the 2x2 matrices which contain the Coulomb wave functions 

All subscripts in the above equations refer to the values of the angular momentum I. Analogously, the 
equation (C.2) has to be replaced by 

x s (r) = F°(r) - mqK s G°(r) , (C.ll) 

where K s is the K-matrix for the strong potential only. 

Matching now the wave function x l ( r ) with x s (r), calculated from the potential V defined in eqs. (C.3) 
and (C.5), at some radius R by equating the corresponding logarithmic derivatives as described above 
one obtains for the K-matrix K l : 

Kl = ^- i F (Fo ~ mqK'Go)- 1 ^ - mqK s G' ) - F'] [G(F - mqK s G G )~ l {F^ - mqK s G' ) - G'] ^ . 

(C.12) 
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The only remaining difficulty is the determination of the matching radius R, because the given solution 
is wrong as long as (C.4) is not valid. On the other hand, it is not possible to extend R to arbitrarily 
large values, because the cosine in eq.(C.5) will cause rapid oscillations. R ~ 10 fm turns out to be a 
good choice, see [7]. We use the value R = 12 fm in the present analysis. 

D Effective range expansion 

In this appendix we collect the formulae for the S-wave effective range expansion. In the simplest case 
of the scattering with the finite-range potential, the quantity fccot(<5o), where 8q is the S-wave phase 
shift and k is the c.m.s. momentum, is well known to have the low-momentum (or effective range) 
expansion: 

fccot^o) = — + \rk 2 + v 2 k 4 + v 3 k 6 + v A k 8 + 0(k 10 ) . (D.l) 

Here a is the scattering length, r the effective range and ^2,3,4 the shape parameters. 

In the presence of the long-range potential the effective range expansion has to be modified. In that 
case one usually defines an effective range function instead of the quantity fccot(<5o), in which the left- 
hand singularities due to the long-range interaction are removed, see reference [50] for more details. 
In the case of the modified Coulomb potential given in eq. (2.53), the effective range function Fq takes 
the form [50] 

F C = C%(v') k cot^) + 2k 7/ h{rf) , (D.2) 



where the quantity rj is given by 
and the functions Cq(t]') (the Sommerfeld factor) and h(rj') read 



v' = ^>, (D.3) 



^V) = e J7_ 1 ' and h w = Re [* (1 + ir]l) ] ~ ln(r?,) ■ (D - 4) 

Here, \£ denotes the digamma function. Notice that the phase shift 5q is the S-wave phase shift of 
the finite-range plus Coulomb potential with respect to Coulomb wave functions. In the notation of 
section 5, 8q should be written as <5^ 1+N . The effective range expansion for the function Fq is 

F c = ~ + -r c k 2 + v$k 4 + v$k 6 + v%k 8 + 0{k w ) . (D.5) 
a z. 

In a general case of an arbitrary long-range interaction, the effective range function may be obtained 
e.g. along the lines of Ref. [100] provided that the long-range potential is weak enough to be treated 
perturbatively. For more discussion on the effective range expansion in presence of electromagnetic 
interaction the reader is referred to [50]. 
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TABLES 



Isospin-symmetric 


Isospin-breaking, finite-range 


Isospin-breaking, long-range 


LO (y = 0): 

static IPE, contact terms 
without derivatives 






NLO (y = 2): 

leading 2PE, contact terms 
with 2 derivatives 


L0(^ = 2): 

M^i / M n o in IPE 


L0(^ = 2): 
static 17-exchange 


NNLO [y = 3): 

subleading 2PE 


NL0(^ = 3): 

isospin breaking in IPE (oc 
eM 2 ), contact term without 
derivatives oc eM 2 




N 3 LO (y = 4): 

sub-subleading 2PE, leading 
3PE, l/m 2 -corrections to 
IPE, 1/m-corrections to 
2PE, contact terms with 4 
derivatives 


NNL0(z/ = 4): 
isospin breaking in IPE (oc 
e 2 /(47r) 2 ), M w ± / M n o in 
2PE, triangle and football 2PE 
diagrams oc C5, ^-exchange, 
m n 7^ m p in 2PE and in the LS 
equation, contact terms with- 
out derivatives oc e 2 /(4-7r) 2 


N 4 L0(z/ = 6): 

l/m 2 -corrections to the static 
17-exchange, 27-exchange. 



Table 1: Dominant contributions to the isospin-symmetric and isospin-breaking parts of the two- 
nucleon force. 
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LEC 


{450, 500} 


{600, 600} 


{450, 700} 


{600, 700} 


°150 
/^np 

WSO 
/inn 
150 

ClSO 

ni 
-^150 

D 2 


-0.0834 
-0.0913 
-0.0880 

1.5007 
-26.9836 

3.7402 


-0.0800 
-0.0892 
-0.0851 
1.8075 
-16.7678 
-2.5565 


-0.1247 
-0.1289 
-0.1272 

2.1217 
-24.7288 

0.8214 


-0.0436 
-0.0544 
-0.0494 
1.8950 
-17.6295 
-2.0771 


C3S1 
C3S1 

-^351 

n 2 

^351 


-0.1498 
0.4144 
-26.3516 
4.8091 


0.1782 
-0.9058 
-13.4902 
2.6661 


-0.1599 
0.6275 
-23.8555 
4.3807 


0.0746 
-0.3557 
-12.4078 
1.8895 



Table 2: The S — wave LECs Cj, Cj and A at N 3 LO for the different cut-off combinations 
{A [MeV], A [MeV]}. The values of the Q are in 10 4 GeV" 2 , of the d in 10 4 GeV" 4 and of 
the A in 10 4 GeV~ 6 . 



LEC 


{450, 500} 


{600, 600} 


{450, 700} 


{600, 700} 


C\p\ 
D 1P1 


0.1862 
2.3257 


0.3374 
1.9180 


0.2072 
2.3968 


0.3444 
1.9213 


C3P0 
-D3P0 


1.1729 
1.0892 


1.2034 
1.2500 


1.1913 
1.2190 


1.2031 
1.4116 


C3P1 
-D3P1 


-0.6334 
4.2369 


-0.6602 
3.8465 


-0.7576 
4.2099 


-0.7193 
3.8756 


C3P2 
D3P2 


-0.5542 
4.1956 


-0.5812 
4.2270 


-0.6217 
4.0340 


-0.6114 
4.1723 


c el 


-0.4516 
2.6303 


-0.2726 
1.7686 


-0.5045 
2.0296 


-0.3352 
1.5516 



Table 3: The LECs Cj and Di in the P-waves and e\ at N 3 LO for the different cut-off combinations 
{A [MeV], A [MeV]}. The values of the d (A) are in 10 4 GeV" 4 (10 4 GeV" 6 ). 
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LEC 


{450, 500} 


{600, 600} 


{450, 700} 


{600, 700} 


D\D2 


-2.2450 


-2.1874 


-2.3398 


-2.2203 


Dzdi 


-1.3988 


-1.7483 


-1.2250 


-1.6620 


D3D2 


-1.4180 


-0.9023 


-1.3578 


-0.8580 


D3D3 


-2.0792 


-1.5493 


-1.7522 


-1.4841 


D e2 


0.2333 


0.2901 


0.2274 


0.2892 



Table 4: The D-wave LECs A at N 3 LO for the different cut-off combinations {A [MeV], A [MeV]}. 
The values of the A are in 10 4 GeV~ 6 . 





Chiral N 3 LO 


Nijmegen PWA 


i 


AZ 


Af 


Aj 


A7 


Af 


A, 




0.359... 0.360 


-0.456... -0.454 


-0.096 ... - 0.094 


0.371 


-0.447 


-0.076 




-0.187... -0.184 


0.191... 0.192 


0.005... 0.008 


-0.186 


0.183 


-0.003 


3 P2 


0.092... 0.093 


-0.031 ... - 0.030 


0.061... 0.062 


0.092 


-0.035 


0.057 


'D 2 


0.014... 0.015 


-0.024... - 0.024 


-0.010... -0.009 


0.014 


-0.023 


-0.009 



Table 5: Effects (in degrees) on the phase shifts 5i, i = { 3 Po, 3 -Pi, 3 P 2 , 1 A}, due to removal of the 
Coulomb interactions (Aj) and subsequently including the pion mass differences in the 1PE potential 
(Af ) at £i ab =_ 10 MeV. The shifts A], Af and A; are_defined as follows: A] = (S pp )i - (S pp )i, 
AJ = (5 n p)j — (<5pp)i> Aj = Aj + AJ = (5 np )i — (S pp )i and (S pp )i denotes the pp phase shifts calculated 
in the absence of the Coulomb interaction. The Nijmegen PWA results are from [38]. The cut-offs A 
and A are varied as specified in eq. (5.1). 
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Chiral N 3 LO 


Nijmegen PWA 


i 


AJ 




Aj 


AJ 


A^ 

i 


Ai 




0.320... 0.325 


-0.789 ... - 0.784 


-0.465 ... - 0.464 


0.342 


-0.785 


-0.443 


3 Pi 


-0.222... - 0.218 


0.293... 0.293 


0.071... 0.075 


-0.221 


0.275 


0.054 


3 P2 


0.185... 0.190 


-0.099... -0.093 


0.088... 0.096 


0.184 


-0.115 


0.069 


'D 2 


0.029... 0.031 


-0.048 ... - 0.048 


-0.020... - 0.017 


0.031 


-0.046 


-0.015 



Table 6: Effects on the phase shifts Si, i = { 3 Pq, 3 -Pi, 3 P 2 , 1 P j 2\, due to removal of the Coulomb 
interactions (Aj) and subsequently including the pion mass differences in the 1PE potential (A^) at 
£i a b = 25 MeV. For notations see Table 5. 





Chiral N 3 LO 


Nijmegen PWA 


i 


AJ 


a: 


A, 


AZ 


Af 


A, 


3 Po 


0.091... 0.109 


-0.879... - 0.866 


-0.775 ... - 0.770 


0.119 


-0.896 


-0.777 


3 Pi 


-0.236... -0.227 


0.320... 0.323 


0.084... 0.095 


-0.233 


0.297 


0.064 


3 P2 


0.251... 0.260 


-0.175... -0.161 


0.078... 0.098 


0.253 


-0.221 


0.032 


l D 2 


0.045 . . . 0.050 


-0.041 ... - 0.041 


0.003... 0.008 


0.049 


-0.034 


0.015 



Table 7: Effects on the phase shifts 5i,i = { 3 Po, 3 Pi, 3 P 2 , 1 P > 2}, due to removal of the Coulomb 
interactions (Aj) and subsequently including the pion mass differences in the 1PE potential (Af ) 
at .Ei a b = 50 MeV. For notations see Table 5. 





NLO 


NNLO 


N 3 LO 


Nijmegen PWA 


a [fm] 


-23.447.. 


. - 23.522 


-23.497.. 


. - 23.689 


-23.585. . 


. - 23.736 


-23.739 


r [fm] 


2.60. 


.2.62 


2.62. 


.2.67 


2.64. 


.2.68 


2.68 


v 2 [fm 3 ] 


-0.46. . 


. - 0.47 


-0.48.. 


. - 0.52 


-0.49.. 


. - 0.51 


-0.48 


f3 [fm 5 ] 


4.3. 


.4.4 


4.0. 


.4.2 


4.0. 


.4.1 


4.0 


t>4 [fm 7 ] 


-20.7. . 


. - 21.0 


-19.9.. 


. - 20.5 


-19.8.. 


. - 20.2 


-20.0 



Table 8: Scattering length and range parameters for the So partial wave using the NLO and 
NNLO potential [21] compared to the N 3 LO results and to the Nijmegen phase shift analysis 
(PWA). The values ^2,3,4 are based on the np Nijm II potential and the values of the scattering 
length and the effective range are from ref. [101]. 
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NLO 


NNLO 


N 3 LO 


Nijmegen PWA 


a [fm] 


5.429. 


. 5.433 


5.424. 


. 5.427 


5.414. 


. 5.420 


5.420 


r [fm] 


1.710. 


. 1.722 


1.727. 


. 1.735 


1.743. 


. 1.746 


1.753 


V2 [fm 3 ] 


0.06. 


.0.07 


0.04. 


.0.05 


0.04. 


.0.05 


0.04 


v% [fm 5 ] 


0.77. 


.0.81 


0.71. 


.0.76 


0.69. 


.0.70 


0.67 


v A [fm 7 ] 


-4.3. . 


. -4.4 


-4.1.. 


. -4.2 


-4.0. . 


. -4.1 


-4.0 



Table 9: Scattering length and range parameters for the 3 5i partial wave using 
the CR NLO and NNLO potential [21] compared to the N 3 LO results and to 
the Nijmegen PWA [93]. 





NLO 


NNLO 


N 3 LO 


Exp 


E d [MeV] 


-2.171... -2.186 


-2.189... - 2.202 


-2.216... -2.223 


-2.224575(9) 


Qd [fm 2 ] 


0.273... 0.275 


0.271... 0.275 


0.264... 0.268 


0.2859(3) 


Vd 


0.0256... 0.0257 


0.0255... 0.0256 


0.0254... 0.0255 


0.0256(4) 




1.973... 1.974 


1.970... 1.972 


1.973... 1.985 


1.9753(11) 


A s [fm- 1 / 2 ] 


0.868... 0.873 


0.874... 0.879 


0.882 . . . 0.883 


0.8846(9) 


P d [%} 


3.46... 4.29 


3.53... 4.93 


2.73... 3.63 





Table 10: Deuteron properties derived from the chiral potential at N 3 LO compared to the NLO 
and NNLO results from [21] and the data. Here, E d is the binding energy, Q d the quadrupole 
moment, i] d the asymptotic D/S ratio, yjr^, the root-mean-square matter radius, As the 
strength of the asymptotic S-wave normalization and P d the D-state probability. The data for 
E d are from [102], for Q d from [103, 104], for rj d from [105] and for As from [104]. For the rms- 
radius we actually show the experimental value for the deuteron "point -nucleon" rms-radius 
from [64]. In the N 3 LO calculation, the cut-offs are varied as specified in eq. (5.1). 
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FIGURES 




r [fm] 



Figure 1: Isoscalar spin-spin 3PE potential using dimensional (DR) and spectral 
function regularization (SFR). The cut-off in the spectral function varies in the 
range A = 500 . . . 700 MeV. 




1 2 3 4 5 6 7 

r [fm] 



Figure 2: The ratio of the isoscalar spin-spin 3PE and 2PE N 3 LO contributions 
using dimensional (DR) and spectral function regularization (SFR). The cut-off in 
the spectral function varies in the range A = 500 . . . 700 MeV. 
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Figure 3: S-wave np phase shifts versus the nucleon laboratory energy. The grid, 
light shaded and dark shaded bands show the NLO, NNLO [21] and N 3 LO results, 
respectively. The cut-offs A and A at N 3 LO are varied as specified in eq. (5.1). The 
filled circles depict the Nijmegen PWA results [38] and the open triangles are the 
results from the Virginia Tech PWA [106]. 
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Figure 4: P-wave np phase shifts and mixing angle e\ versus the nucleon laboratory 
energy. For notation see Fig. 3. 
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50 100 150 200 250 



Lab. Energy [MeV] 

Figure 5: 3 Pi-wave np phase shift versus the nucleon laboratory energy. The solid 
and dashed lines correspond to the LECs C%pi and D^px from the first and second 
lines in eq. (6.3), respectively. For remaining notations see Fig. 3. 
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Figure 6: D-wave np phase shifts and mixing angle e% versus the nucleon laboratory 
energy. For notation see Fig. 3. 
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Figure 7: F-wave np phase shifts and mixing angle €3 versus the nucleon laboratory 
energy. For notation see Fig. 3. 
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Lab. Energy [MeV] 

Figure 8: G-wave np phase shifts and mixing angle versus the nucleon laboratory 
energy. For notation see Fig. 3. 
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e [deg] 



Figure 9: np differential cross section and vector analyzing power at £j a b = 25 MeV. 
The Nijmegen PWA result is taken from [107]. Data for the cross section are taken 
from [108] and for the analyzing power from [109, 110]. The cut-offs A and A are 
varied as specified in eq. (5.1). 
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e [deg] 



Figure 10: np differential cross section and vector analyzing power at £i a b = 50 
MeV. Data for the cross section are taken from [108, 111] and for the analyzing 
power from [112, 113, 114, 115, 110]. For remaining notations see Fig. 9. 
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Figure 11: np differential cross section and vector analyzing power at E\ a \, = 96 MeV. 
Data for the cross section are taken from [118, 119, 116]. Data for the analyzing 
power are at £?i a b = 95 MeV and taken from [117]. For remaining notations see 
Fig. 9. 
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Figure 12: np differential cross section and vector analyzing power at E\ a \y = 143 
MeV. Data for the cross section are at = 142.8 MeV and taken from [120] and 
for the analyzing power from [121]. For remaining notations see Fig. 9. 



61 



0.7 



NLO 

NNLO 

NNNLO 




0.15 - 



M 01 







1 1 1 


i 1 i < i < 






'I s\ 




/> / \ 




■if \ 




ii \ 




i' / \ 




■/ / ^ 




// / 




/ 




/ 




// / 




''/ 






/ 




/// 








" V 




^ . i 


1,1,1, 







r [fm] 



8 



10 



Figure 13: Coordinate space representation of the S- (upper panel) and D-wave 
(lower panel) deuteron wave functions at NLO, NNLO and N 3 LO for the cut-offs: 
A = 550 MeV, A = 600 MeV. 
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